{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = np.genfromtxt(\"Troe.dat\", dtype=str)\n",
    "Header = (data[0,:]).tolist()\n",
    "solTchem = (data[1:,:]).astype(np.float)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "campsol = np.loadtxt(\"troe_results.txt\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "A_indx = Header.index('A')\n",
    "B_indx = Header.index('B')\n",
    "C_indx = Header.index('C')\n",
    "t_indx = Header.index('t')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(46, 9)\n"
     ]
    }
   ],
   "source": [
    "print(np.shape(data))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x11d0675d0>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3dd3wUdfrA8c+zLYUESCChBQi9HVIVERHEcoAFu3B6iHr2epY729l/d3rqeZyn52GvoKIoKggowomigkgTKaEHREJIgEDa7n5/f8wmhJBkN8lONsk+79drXrszOzvzLBvm2W+Z71eMMSillIpejkgHoJRSKrI0ESilVJTTRKCUUlFOE4FSSkU5TQRKKRXlXJEOoLpatmxp0tPTIx2GUko1KD/88MMeY0xKRa81uESQnp7O0qVLIx2GUko1KCKytbLXtGpIKaWinCYCpZSKcpoIlFIqymkiUEqpKKeJQCmlopxtiUBEXhaR3SKyupLXRUT+JSIZIrJSRAbaFYtSSqnK2VkieBUYXcXrY4BugeVq4D82xqKUUqoStiUCY8z/gL1V7DIOeN1YvgWai0gbu+JZsmUvT81dR7HPb9cplFKqQYrkDWXtgO1l1jMD234pv6OIXI1VaqBDhw41OtmyrTk8Mz+Da0Z0we3UphGlVO35/D6K/EUU+4sp8hVR7Cs+/NxfbL0W2Ob1e63FeA8/93tLX6tqH5/x4fV7GdtpLINbDw7754hkIpAKtlU4S44xZgowBWDw4ME1mknHFbj4+/w6EY9SjZExhgJfAfnefA4VH7IevdZjfnE++d58Cn2FpUuBt6Ba6yUX99KLvb8Iv7GvhsElLlwOa3E6nDjFSb+UfgymcSWCTKB9mfU0YKddJ3M5rLzj1aohpeoNYwz53nwOFB0grziv9DGvKI8DxQesxzLbDhYfPHxxL3PBL1lMxb8lK+UUJzHOGGJdsXicHmKdscQ4Y6zFFUOSO4kYZwwep4cYZwxuhxuP04PH4cHlcFnPnR5ru8OD2+ku3afso9vhLr2olyxuOXpb6SIuRCr6rWyPSCaCmcCNIjINGALsM8YcVS0ULi6n9Y+qJQKlws8YQ15xHrkFueQW5pJTmENuYe4R6/sK95FTkMP+ov3kFeWRV2xd2H3GV+WxHeKgibsJie5E4t3x1uKKJzk2mXh3PHGuOOJcccS7rOcl20rWS7bFumJLL/SxLuvR5Whww63ZwrZ/BRGZCowEWopIJvAA4AYwxjwPzALGAhnAIeByu2KBwyWCYk0ESoXE5/eRU5jD7kO7yTqURVZ+FlmHstidv5vcgiMv7vsK9+E13gqP4xQnzWKakRSTRLOYZrRNaEuiO5EETwIJ7gQSPdbzstsS3AnWNk8i8a74Ov11HI1sSwTGmAlBXjfADXadvzyXI9BG4NNEoKKbMYa9BXvJys+q8CJfsp6dn13hr/Xk2GSSYpJoHtuc9Kbp9EvpR1JsEs1jmtM8pvkRz5vHNifBnYBDtINGfRY15aKSqqFiv7YRqMYv35vPjgM72JG3g8y8TDIPBJa8THbk7SDfm3/Ue5JikmgZ35LUuFS6JXUjJS6F1PhUUuJSSIm3nreIbYHb6Y7AJ1J2ip5E4NBeQ6rx8Bs/uw/tPuLiXvZivyd/zxH7x7niSEtMo31ie4a2HUq7hHalF/nU+FRaxrXE4/RE6NOoSIuaROAsaSPQXkOqgcnOz2ZD7gYycjLIyM1gQ84GMnIzOOQ9VLqPQxy0im9FWmIaw9sNJy0xjbSENNIS02iX0I7k2GStZ1eVippE4NZeQ6qeO1h88IgLfUZOBhtyN7C34PAN+s1jmtMtqRvndD2Hzs060z6xPWmJabRp0karbFSNRU0iOFwi0ESgIu/Xg7+yPGs5P2f/XHrx33nw8G00ca44ujbvyoi0EXRL6kbX5l3pltSNFrEt9Je9CruoSQQlw0roDWWqrnn9XtbnrOfH3T+yYvcKlmct55eD1i0zLnGR3szqeXN+9/Pp1rwbXZO60i6hnfa0UXUmahJBSYnAq1VDymb7CvexMmsly7OWs3z3clbtWVXaSyc1PpUBqQOY2Hsi/VP70yOph1bpqIiLmkRQ0kagiUCFkzGGbQe28ePuH1m+ezkrslaQkZsBWDdSdU/qzjldz2FA6gD6p/SndZPWWrWj6p2oSQQl3Ue1akjV1qHiQyzeuZiFmQv5X+b/yC7IBiDRk0i/lH6MTh9N/9T+9G3Zl3h3fISjVSq46EkETm0sVjW3M28nCzMXsnD7Qr7f9T3F/mIS3YkMazeMIW2G0D+lP52bd9Z6fdUgRU0iKG0s1juLVQh8fh+r9qyyLv6ZC9mQswGAjk07MqHnBEakjWBAqwG4HVq/rxq+qEkEh4eh1hKBqlheUR6Lf1nMgu0LWLRjEXsL9uIUJwNbDeSOwXcwIm0E6c3SIx2mUmEXNYmgpESgdxarsvYW7GX25tks3L6QJb8uwev30tTTlBPbncjI9iM5oe0JNItpFukwlbJV1CQCl/YaUgF+4+fbX77l/fXvM3/7fLx+L52adeL3vX7PSWkn0T+1v45Tr6JK1Py1a68htevgLj7M+JAPMz5kR94Omsc0Z0LPCZzX9Ty6JnWNdHhKRUzUJAK39hqKSsX+Yr7K/Ir3N7zPoh2L8Bs/x7c5nlsH3sqoDqN0xE2liKJE4NJeQ1Fl2/5tfLDhAz7a+BF78veQEpfClb+5knO7nUv7xPbBD6BUFImeRKCDzjV6hb5CPt/6OR9s+IDvd32PQxyc1O4kzu9+Pie2O1Hr/ZWqRNT8zyhJBDoMdeOTeSCTN39+k483fsz+ov20S2jHTQNuYlyXcbRq0irS4SlV70VNIigddE4bixuNXQd3MWXlFGZsmIGIcEqHUziv23kMaTNE7/BVqhqiJhGICG6nUKwlggZvT/4eXlr1Eu+uexc/fi7ofgFXHXMVqfGpkQ5NqQYpahIBWF1ItUTQcOUW5PLyTy8zbe00inxFjOs6jmuOuYa2CW0jHZpSDVp0JQKnaGNxA3Sg6ACvr3mdN9a8waHiQ5zR+Qyu63cdHZp2iHRoSjUKUZUI3E6Hdh9tQA4VH+LttW/zyupX2F+0n9M6nsb1/a7Xm7+UCrMoSwSig841AAXeAt5d9y4vrX6JvQV7GZE2ghv630CvFr0iHZpSjVJUJQKXw0GRthHUW8W+Yt7f8D4vrHyB3fm7Ob7N8dw44Eb6pfSLdGhKNWpRlQg8LkfwEsEvK2Hh4zBoEnQ7rU7iinbGGD7b8hn//OGf7Dy4k4GpA3nspMc4tvWxkQ5NqagQVYnA5ZCqh6H+4VWYdSf4imHtJzD8dhj1F9A5Zm2TnZ/No98+yufbPqdXci8eGPoAQ9sO1Xl9lapDUZUI3E5H5Ymg8ICVBNoPgXOfhwV/g6+egpY9oN/FdRtolJi3dR6PLH6EvOI8bh14K5P6TMLpcEY6LKWiTnQlApej8u6jGZ+DrwhOvgeapcFZz8DutTDnbquKKD65boNtxPYV7uP/vvs/Zm+eTa/kXvzfif9Ht6RukQ5LqagVVffhu6uqGlo3G+KSIe04a93hgLMmQ34uzLu/7oJs5BZuX8g5H53DvC3zuL7/9bx1xluaBJSKMFsTgYiMFpF1IpIhIndV8HoHEflSRH4UkZUiMtbOeNzOShqLfV5YPwe6jwZnmUJS69/AkGth+VuQvdHO0Bq9A0UHuG/Rfdw4/0aSYpN4+4y3ua7fdTr5u1L1gG2JQEScwLPAGKA3MEFEepfb7T7gXWPMAGA88Jxd8YB1Z3GF3Ue3LYaCXOgx5ujXht0CTg8s+oedoTVq3+z4hnM/OpePN33MVX2vYtoZ0/SeAKXqETtLBMcBGcaYTcaYImAaMK7cPgZoGnjeDNhpYzx4KmssXv+ZdbHvMuro1xJbwcDLYMU0yN1mZ3iNzsHigzy8+GGu+fwa4t3xvDnmTW4eeLPOCqZUPWNnImgHbC+znhnYVtaDwKUikgnMAm6q6EAicrWILBWRpVlZWTUOyFXZncW/rIA2/SAmoeI3DrsZEPjmmRqfO9os2bWE82eez/T107ms92W8e+a79E3pG+mwlFIVsDMRVNQRvPxVeALwqjEmDRgLvCFy9EDyxpgpxpjBxpjBKSkpNQ6o0u6jOVshqVPlb2yWBn0vhB/fgoJ9NT5/NMj35vP4949zxZwrcIiDV0e/yh3H3kGsKzbSoSmlKlFl91ERuS2EYxw0xvy3gu2ZQNnJYdM4uurnSmA0gDFmsYjEAi2B3SGct9o8TgfF5Qed8xXD/kxISq/6zUOugRVvw49vwtAb7AivwVuRtYJ7F93L1v1bmdBzArcOvJV4d3ykw1JKBRGsRHAnkAAkVrHcXsl7lwDdRKSTiHiwGoNnlttnG3AKgIj0AmKBmtf9BOFyCsXecoWSfdvB+IMngrb9ocNQ+O6/4PfZFWKDNWvTLCZ9NokiXxEvnP4C9wy5R5OAUg1EsBvK3jDGPFzVDiLSpKLtxhiviNwIzAGcwMvGmJ9E5GFgqTFmJlYSeUFE/ohVbTTJGGPb8KAVDkOds8V6TOoY/ABDroH3JlldTXva2tO1wTDG8NpPr/HUD08xqNUgJp88mWYxzSIdllKqGqpMBMaYPwU7QFX7GGNmYTUCl912f5nna4BhwcMMD7fTQZG3skSQHvwAPc+CxLaw5EVNBIDP7+OJpU/w1s9vcXrH0/nr8L8S44yJdFhKqWqqsmpILBeJyIWB56eIyL9E5PqKGnXrO3dFM5TlbLG6jia2CX4Ap8salXTjF1F/g1mhr5A7/3cnb/38Fpf2upQnRjyhSUCpBirYxfxZ4CLg98AbwLXAUuAk4Gl7Qwu/iquGtkLzDhDqYGcDJ4LDBT+8Ev4AG4h9hfu4eu7VzNs6jzsG38Gfj/szjob3u0ApFRCsjWC4MaaviLiBXUAbY0yRiLwN/Gh/eOHlclqDzhljDg9znLMFmofQPlCiaRvoeYbVe+jke8EdZ0us9dXOvJ1c9/l1bD+wnSdOeoLRnUZHOiSlVC0F+xnnBTDGFANLAncIY4zxAg2u64zHaV38vf4y1UM5W0JrHyjr2D9Afg6s+ShssTUE6/au49JZl5J1KIv/nvZfTQJKNRLBEsEuEUkAMMaU/q8XkdZAkZ2B2cHltD5u6U1l+bnWGEPVTQTpw6FFV1j6cngDrMcW71zMZZ9dhkMcvDbmNZ09TKlGpMpEYIwZY4zJq+Cl/cCZ9oRkH3dpIgiUCHK3Wo/VTQQiMPgK2P4d7FodvgDrqY83fsz1n19P24S2vDn2TR02WqlGJuQWPhE5RkTOFpHzsO4GPtG+sOzhDlQNlZYIqnMPQXn9JoAzplE3GhtjeHHVi9yz6B4GthrIa6Nfo3WT1pEOSykVZiHNUCYiLwPHAD8BJd1uDPCBTXHZwl2+amjfDuuxWftK3lGF+GT4zXmw4h049UGISQxLjPWFz+/jse8fY9q6aYxJH8OjJz6qo4Yq1UiFOlXl8caY8nMJNDiliaBkmIn8HEAgtnnNDjj4SlgxFVa9Z1UVNRIF3gLu+uouvtj2BZf3uZxbB92q3UOVasRC/d+9uIJJZRocj8v6uEW+QIenglyIbWpNS1kTaYOhVV9Y8jLYNzJGncotyOWquVcxf9t87jruLm4bfJsmAaUauVD/h7+GlQzWBaaUXCUiK+0MzA6eQImg0Fum11BNSwNgNRofewX8ugoyl4YhwsjaW7CXiZ9NZE32Gp4a+RSX9Lok0iEppepAqFVDL2PdXbyKw20EDU5MSYmgJBEU5EJcUu0O2vdCmHu/Nf5Q+4bbpbLAW8DN829mZ95Onj/tee0eqlQUCbVEsM0YM9MYs9kYs7VksTUyG3jKJ4L8XIirRYkArEbifhfDTzPgYHYtI4wMv/Fz39f3sSJrBX8b/jdNAkpFmVATwVoReVtEJojIeSWLrZHZ4HAbQZkSQW2qhkoc+wfwFcLyN2t/rAj494//Zs6WOdw26DZO63hapMNRStWxUBNBHFAInA6cFVga7A1lh0sEObUvEQCk9oKOJ8KSl6D8oHb13IwNM3hh1Qtc0P0CJvWZFOlwlFIREGyqygnAXGPM5XUUj608Ze8jMKb2jcVlHXslTL8cMj6H7qeH55g2+/aXb3l48cOc0PYE7hlyz+GB+JRSUSVYiaAj8J6IfCUiD4rIEGnAV4uSqqFCrx+KD4G/ODwlAoBeZ0FCa/h+SniOZ7ONuRu57cvbSG+WzpMjnsTtcEc6JKVUhAQba+gxY8woYCywArgCWBZoL5goIq3qIshwOaLXUH6utTFcJQKnGwZfDhnz6v2kNXvy93DDFzfgcXp49pRnSfQ0rruilVLVE1IbgTHmgDFmhjHmGmPMAOBRIAV43dbowuyIxuKCQCIIV4kAYNDl4HDD9y+E75hhVuAt4Jb5t5Cdn82/T/k3bRPaRjokpVSEBWsjGFjFy18aY54Kczy28jgrKBHU9j6CshJbQe9xsPwtGHUfxCSE79hh4Dd+7ll0D6v2rOLpkU/zm5a/iXRISql6INgNZVVd6A0wKoyx2O6I+wjyc6yN4aoaKjHkGlg93RqD6LirwnvsWpq8bHLp9JKndDwl0uEopeqJKhOBMebkugqkLhzRfdSOqiGAtGOh3SD47nlrULqajmMUZtPXT+fl1S9zUfeLmNh7YqTDUUrVIyFdpUTELSI3i8j0wHJjYB7jBuWI+QjC3VhcQgSGXAfZGbDxi/Aeu4a+2fkNj377KMPaDePuIXdrN1Gl1BFC/bn6H2AQ8FxgGRTY1qCICB6Xg8LSxmKBmKbhP1HvcVZX0m+fC/+xqykjJ4PbF9xO5+adefKkJ3E5Qh1eSikVLUK9KhxrjOlXZn2+iKywIyC7xTgdVtWQPxdim9lTdePywHF/gPmPwu6frTuPI6Ckm2icK47nTnmOBE/9arxWStUPoV4FfSLSpWRFRDoDPntCspfH5TjcRhDu9oGyBl8JrjhY/Kx956hCvjefm764iZzCHJ455RmdYlIpValQE8GdwJciskBEFgLzgdvtC8s+pYkgPye8XUfLi0+G/r+Dle9C3m77zlMBv/Fzz1f38FP2Tzw2/DH6tOhTp+dXSjUsod5Q9gXQDbg5sPQwxnxpZ2B2cTsd1g1l4RxnqDLHXw++Imuugjr0zx/+yefbPufOY+9kVIcG1cNXKRUBofYacgK/BUYCpwA3iMhtNsZlG4/LYfUasrtqCKBlV+gx1rrTuOigvecKmL15Nq/89Arje4zn0l6X1sk5lVINW6iNxR8DBTTwGcrAuru49M5iu0sEAMNuhnWfwo9vwZCrbT1V1qEsHv32UY5JOYY/H/dn7SaqVBgVFxeTmZlJQUFBpEOpUmxsLGlpabjdoffwDzURpBljjqluQCIyGpgMOIEXjTGPVbDPRcCDWHcqrzDG/K6656kOj8tBYbGvbkoEAB2Oh/ZDYPEzMPgKcNrTfdMYw0OLH6LQV8ijwx7VbqJKhVlmZiaJiYmkp6fX2x9Zxhiys7PJzMykU6dOIb8v1Mbi2SJSrUH2A9VJzwJjgN7ABBHpXW6fbsDdwDBjTB/g1uqcoyY8LgdSfAj83ropEQAMuwVyt8GaD207xUcbP2Jh5kJuGXgLnZqF/geglApNQUEBLVq0qLdJAKx7pVq0aFHtUkuoieBbYIaI5IvIfhE5ICL7g7znOCDDGLPJGFMETAPGldvnKuBZY0wOgDHG9u41MS4HMd591kpdlAgAuo+Blj1g0dPWhDhhtuvgLh7//nEGpg7kkl6XhP34SilLfU4CJWoSY6iJ4ClgKBBvjGlqjEk0xgS7JbcdsL3MemZgW1ndge4i8rWIfBuoSjqKiFwtIktFZGlWVlaIIVfM43TgKT5grcQ2q9WxQuZwwIl/hF9Xw4a5YT20MYYHvnkAn/Hx6LBHcUj9GNtIKdVwhHrV2ACsNqZaP2crSkvl3+/C6pY6EpgAvCgiR/1MN8ZMMcYMNsYMTklJqUYIR3M7HTi9+dZKXd5p2/cCaNYBvnoqrKWC99a/xzc7v+G2QbfRvmn7sB1XKVU/zZgxAxFh7dq1YTtmqIngF2CBiNwtIreVLEHekwmUvTKlATsr2OcjY0yxMWYzsA4rMdjG43Lg9Afqz9xxdp7qSE631YNo+3ew9euwHDLzQCZPLn2SIW2GcFGPi8JyTKVU/TZ16lROPPFEpk2bFrZjhpoINgNfAB4gscxSlSVANxHpJCIeYDwws9w+HwInA4hIS6yqok0hxlQjHpcDpzcCiQBgwKXQJBUW/r3Wh/IbP/d/cz8OcfDICY9olZBSUSAvL4+vv/6al156KayJIKQ+hsaYh6p7YGOMV0RuBOZgdR992Rjzk4g8DCw1xswMvHa6iKzBGrvoTmNMdnXPVR0elwN/aYmgiZ2nOpo7Dk68FebcA1u/gY4n1PhQU9dOZcmuJTx0wkO0SWgTxiCVUsE89PFPrNkZrL9M9fRu25QHzqp6OJgPP/yQ0aNH0717d5KTk1m2bBkDB1Y1kWRoqvwZKSIPBjtAVfsYY2YZY7obY7oYY/4vsO3+QBLAWG4zxvQ2xvQ1xoQvxVXC43Tg9gfaCOq6RADWvMZNUmHBUbdUhGzr/q3884d/MrzdcM7tem4Yg1NK1WdTp05l/PjxAIwfP56pU6eG5bjBSgR/CNJNVLCqfB4MSzR1IMbloMhXYJVR3PF1H4An3rqvYO69sHUxdBxarbf7/D7uW3QfbqebB4Y+0CC6synV2AT75W6H7Oxs5s+fz+rVqxERfD4fIsLf//73Wl8HglUsv8CRbQLll4TAPg2Gx+XAYyLURlBi8BXQJAUWVr9U8MaaN1ietZy7j7ubVk1a2RCcUqo+mj59OhMnTmTr1q1s2bKF7du306lTJxYtWlTrYwebs7jabQP1XazbicMUWSuRSgSlpYL7qlUq2Ji7kWd+fIaT25/MmZ3PtDlIpVR9MnXqVO66664jtp1//vm8/fbbDB8+vFbHjroBaWJcDpxShHHGIA5n5AIZfAV8PRkW/A0uK9+Z6mhev5f7Ft1HvDue+4fer1VCSkWZBQsWHLXt5ptvDsuxo67PYZzHSSyFGFeESgMlPE1g2K2weSFsWhh095dXv8zq7NXce/y9tIxrWQcBKqWiRdQlgliXkziK8Ec6EQAceyU0bQdfPFTl3cbr9q7jPyv+w2/Tf8vo9ApH4VBKqRoLqWpIRFKwBohLL/seY8wV9oRlnziPE68U4nfFRjoUq41i5F0w8yZY+yn0Orrev9hfzH1f30dTT1PuHXJvBIJUSjV2oZYIPgKaAZ8Dn5ZZGpxYt4M4ivA660GJAKDf76BFN/jiYfD7jnr5nbXvsHbvWv5y/F9IirVxjmWlVNQKtbE43hjzZ1sjqSOxLic+CvHVl0TgdMGo++C9y2DFNBhweBjpnIIcnlvxHEPbDOWUDqdEMEilVGMWaongExEZa2skdSTW4yReCil2xEQ6lMN6j4M2/a0eRN7C0s3PLn+WQ8WH+NOxf9JeQkop24SaCG7BSgYFgUlpQpmYpl4qaSwudtSDNoISInDqg7BvuzXRPbA+Zz3vrX+Pi3pcRNekrhENTylVPzidTvr370+/fv0YOHAg33zzTViOG+qgc8FGGm0w4jxODIUU1acSAUCXk6HrqbDw75i+F/H494+T6Enkhv43RDoypVQ9ERcXx/LlywGYM2cOd999NwsXBu9+HkzI3UdF5GwReTKwNNjbWmPdDuKkiCKpRyWCEr/9KxQfZP7cW/l+1/fc0P8GmsXU0SxqSqkGZf/+/SQlhacDSajdRx8DjgXeCmy6RURONMbcVcXb6qVYlxOhkBzqWYkAIKUHhcdeyZM7PqZr83Qu7H5hpCNSSlVk9l2wa1V4j9m6L4ypevyx/Px8+vfvT0FBAb/88gvz588Py6lD7TU0FuhvjPEDiMhrwI9Ag0sEcR4nQhEFUg8TAfBGq/Zk7nYx5UARLongEBhKqXqnbNXQ4sWLmThxYulopLVRnbGGmgN7A88bbH1FjMMgUkxBPSwRZB3K4oWf3+TkxM4MXbkA1nwEfc6JdFhKqfKC/HKvC0OHDmXPnj1kZWWRmppaq2OF2kbwN+BHEXk1UBr4Afhrrc4cIRKYpjK/HiaCycsmU+Qv4o5RT0NqH5j3FyjOj3RYSql6aO3atfh8Plq0aFHrY4Xaa2iqiCzAaicQ4M/GmF21PnskBC6sh4wnwoEcae3etXy08SMu73M5HZp3htF/g9fPhsX/hpPujHR4Sql6oKSNAMAYw2uvvYbTWfsq5CoTgYj0NMasFZGSSTEzA49tRaStMWZZrSOoa8UHATjkr1+JYPKyyTT1NOUPx/zB2tB5BPQ6C776B/S/BJq2jWyASqmI8/mOHoYmHIKVCG4DrgaequA1A4wKe0R2C5QIDtajEsHSXUtZtGMRtw26jaaepodfOO0RWD8HPn8Izvtv5AJUSjVqwWYouzrwdIwxJfM7WkTqY0f8EBQfAiDP745wIBZjDJOXTSY1LpUJPScc+WJyJxh6Iyz6hzWRTYchkQlSKdWohdpYXNF9zOG5t7mulZQI6kkiWJi5kOVZy7m2/7XEVjQ09vDboGkafHSDNhwrpWxRZSIQkdYiMgiIE5EBIjIwsIwE4uskwnALXEwP+CKfCHx+H5OXTaZj046c07WSbqIxiTDu35C9AeY/WrcBKqWiQrA2gt8Ck4A04B9lth8A7rEpJnsFqob214NEMGvzLDJyM3jipCdwO6qIp8vJMPhKWPws9DwDOp5Qd0EqpRq9YG0ErwGvicj5xpj36ygmewVKBPt8kW0s9vq9PL/ieXom9+T09NODv+G0h2HjF/Dh9XDd19acx0opFQYhtREYY94XkTNE5E8icn/JYndwtgiUCPYVV+em6vCbtXkW2w5s49p+1+KQEL6GmAQY9xzkbIZ5D9gfoFKqXtq1axfjx4+nS5cu9O7dm7Fjx7J+/fpaHTOkRCAizwMXAzdh3VB2IdCxVmeOlJISgTdyicDr9zJl5RR6JvdkVPtq9MBNHwZDroMlL8CmBbbFp5Sqn4wxnHvuuYwcOZKNGzeyZs0a/vrXv/Lrr7/W6rih9ho6wRgzEcgxxjwEDAXa1yI4JmsAACAASURBVOrMkRIoEeRGsEQwe/Nstu7fyrXHXFv9waJOuR+Su8BHN0JBg5wbSClVQ19++SVut5trr722dFv//v0ZPnx4rY4b6tWw5B6CQyLSFsgGOtXqzJFSnI8fJ3nekKdiCKuS0kCPpB6c3OHk6h/AEw/nPg8v/xbm3gtnPxP+IJVSVXr8+8dZu3dtWI/ZM7knfz6u6qnhV69ezaBBg8J6Xgi9RPCxiDQHngCWAVuAqWGPpi54C/E53Pj8hmKfv85PP3vzbLbs3xJ620BF2h8HJ9wEy16HDfPCG6BSKuoELRGIiAP4whiTC7wvIp8AscaYfbZHZwdfET6H1WMov9iH21l3JQOf38eUlVPontSdUR1qOTrHyHus4Sdm3gTXL4a48MxUpJQKLtgvd7v06dOH6dOnh/24Qa+CgcloniqzXhhqEhCR0SKyTkQyRKTSSWxE5AIRMSIyOKSoa8NbiN9pJYKCYnsGcKrMvG3z2LJ/C9ccc03NSwMl3LFWFVHebvjs7vAEqJSq10aNGkVhYSEvvPBC6bYlS5bUet7iUK9Gc0XkfKlGy6aIOIFngTFAb2CCiPSuYL9E4Gbgu1CPXSu+IkygRFBQVHdVQ8YYXlz5IulN0zm146nhOWjbATD8dlgxFdbOCs8xlVL1logwY8YM5s2bR5cuXejTpw8PPvggbdvWbnTiUBuLbwOaAF4RKcDqQmqMMU2reM9xQIYxZlPgA0wDxgFryu33CPB34I7qBF5j3gJMSYnAW3clgkU7FrEuZx2PDHuk9qWBsk66E9bPho9vgQ7HQ3xy+I6tlKp32rZty7vvvhvWY4Z6Q1miMcZhjPEYY5oG1qtKAgDtgO1l1jMD20qJyACgvTHmk6oOJCJXi8hSEVmalZUVSsiV8xZhnNbsZPlFdZcIXlz1Iq2btOaMTmeE98AuD5zzH8jPgU9vD++xlVJRIdQbyr4IZVv5XSrYZsq83wE8DQS9ehljphhjBhtjBqekpATbvWq+QqjjNoLlu5ezbPcyJvWZhNtpwxhHrfvCiD/DTx/ATzPCf3ylVKMWbPTRWBFJBlqKSJKIJAeWdCBYpVQmR950lgbsLLOeCPwGWCAiW4DjgZm2Nxh7C8FtlQgKvHXTRvDK6ldoFtOMc7uea99JTvyj1WbwyW1WA7JSKuyMMcF3irCaxBisRHAN1kT1PQOPJctHWA3BVVkCdBORTiLiAcYDM8sEu88Y09IYk26MSQe+Bc42xiyt9qeoDl8R4iqpGvLaeiqATfs28eX2LxnfYzzxbhtH7na64JznoeggfPJHaAB/sEo1JLGxsWRnZ9frZGCMITs7m9jY6s0bFmz00cnAZBG5yRhTrVtYjTFeEbkRmAM4gZeNMT+JyMPAUmPMzKqPYBNvIQ5PAgB5hfZXDb3+0+t4nJ6jZx+zQ2pPGHUvzLvfGrL6hBvtP6dSUSItLY3MzExq3U5ps9jYWNLS0qr1npB6DRljnhGRE4D0su8xxrwe5H2zgFnltlU4aqkxZmQosdSarwiXx8qWeQXFtp5qT/4ePt74MeO6jqNFXAtbz1Vq6E2QuRTm3gdJ6dDrzLo5r1KNnNvtplOnhjmyTjChNha/ATwJnAgcG1jsv/nLDt5CnIE2ggMF9lYNvbvuXYr8Rfy+9+9tPc8RHA44bwq0GwgfXAU7ltXduZVSDVKo9xEMBnqb+lw5FipfEU53LDEuBwcK7UsEhb5C3ln3DiPSRtCpWR3/inDHwYRp8MIpMHU8/OELaN4wB4tVStkv1DubVgOt7QykzngLweUhMdbNARurhj7d9Cl7C/YysfdE285RpYRUuOQ9KC6Aty/SIauVUpUKNRG0BNaIyBwRmVmy2BmYbXyF4IyhaayL/TZVDRljeGPNG/RI6sGxrY+15RwhSe0JF70Ge9bDe5eBz942EaVUwxRq1dCDdgZRp7xFgRKBizybEsGSXUvIyM3g4RMerv7EM+HW5WQ482lrlNJZd1rPIx2TUqpeCbXX0EIR6Qh0M8Z8LiLxWF1CGxZjwFsAzhhbq4be+vktkmKSGNt5rC3Hr7aBE2HvJlj0NLToYs1loJRSAaH2GroKmA78N7CpHfChXUHZxu8FDLhiSIhx2dJrKPNAJgsyF3BB9wuICYxpVC+Muh96nwNz/wJrGmatnlLKHqG2EdwADAP2AxhjNgCpdgVlG2+h9ei0qobsSATvrnsXQbiox0VhP3atOBzW/AVpg+GDqyHzh0hHpJSqJ0JNBIXGmKKSFRFxUWYAuQbDF/gIrlhbqoYKvAV8kPEBozqMonWTetjJyh0H46daPYqmXgw5WyMdkVKqHgg1ESwUkXuAOBE5DXgP+Ni+sGxSUiIINBYfLPLh84cvn3225TP2Fe5jfI/xYTtm2CWkWN1KvUWBbqUNc8ZRpVT4hJoI7gKygFVYA9HNAu6zKyjb+EqqhmJIjLXayfPCeFPZtLXT6NKsS2S7jIYipQdc/AZkZ8C72q1UqWgXaiKIwxo07kJjzAXAy4FtDYu3pGrocCIIV/XQT3t+4qfsn7i458WR7zIais4j4KzJsOlL+PQ2Ha1UqSgWaiL4giMv/HHA5+EPx2a+so3F1gQx4Wowfnf9u8S54jir81lhOV6dGHApDL8Dlr0OX0+OdDRKqQgJ9YayWGNMXsmKMSYvcC9Bw1K2ROAuKRHUPhHsL9rP7M2zGdtpLAmBIa4bjJPvte4x+PwBa7TSPudEOiKlVB0LtURwUEQGlqyIyCAg356QbFRhiaD2VUOfbPyEfG8+F/a4sNbHqnMOhzXncfshMOMa2L4k0hEppepYqIngVuA9EflKRL4C3gEa3qwnpb2GwtdYbIzhvfXv0btFb/q06FPbCCPDHQvj34bE1tZopdkbIx2RUqoOhZQIjDFLsKarvA64HuhljGl4dySV3EdQptdQbQeeW7lnJRm5GVzQ/YLaRhdZTVrCJdPB+OGl02FHw/t6lVI1E2qJAKzJaI4BBgATRCRC4yvXQtn7CGLCUzX0wYYPiHPFMbZTPRlXqDZadoMr54InHl49E9bPjXRESqk6EF0zlHkP30cQ63bgckitGosPFh9k9ubZjE4fTRN3kzAFGWEtu8GVn1uPU8dbPYqUUo1adM1Q5jtcIhCRwHhDNS8RzN48m3xvPud3Pz9MAdYTia1g0qfWzWYzb4J9O2DkXTp8tVKNVHTNUFamjQAgMdZdqzkJZmyYQZdmXTim5THhiK5+iUmE370D/X4HCx+zEoLP3jmelVKREWqJoGSGsu+BwpKNxpizbYnKLiUXMqcHoFYjkG7M3cjKPSu5Y/AdDeNO4ppwuuGc56BZO/jfE5D3K1zwCsQ0sHsllFJViq4ZykpLBFZDcW3mJPgw40Nc4uLMzmeGK7r6SQRG3QdN28Knt8NrZ8Lv3rMGr1NKNQqhdh9dCKwFEgPLz4FtDUu5RJAY62Z/DdoIiv3FzNw4kxHtR9AirkU4I6y/Bl9h3Wuwey28dKrea6BUIxJqr6GLgO+BC4GLgO9EpOF1nPcfWTXUtIZVQ4syF7G3YC/ndI2y4Rh6jIFJn0DhAXjpNMhcGumIlFJhEGpj8b3AscaYy4wxE4HjgL/YF5ZNfEUgDnBY0y0nxrpqdGfxRxs/Ijk2mWHthoU7wvovbTBcOc9qTH71TFg3O9IRKaVqKdRE4DDG7C6znl2N99YfviJwuEtXE2Pd5BV6qU6v2NyCXBZmLuSMzmfgLnOsqNKii5UMUnvCtN/B0pcjHZFSqhZCvZh/JiJzRGSSiEwCPgUa3k9Bn7e0WgggIdaFz284VOQL+RCzt8zG6/cyrss4OyJsOBJS4bJPoOup8MkfYf6jOqeBUg1UqI3FdwL/xRpioh8wxRjzJzsDs4WvqLShGCAp3nq+92BRZe84ysyMmfRI6kGP5B5hD6/BiUmw5kAe8Hure+lHN+hsZ0o1QFUmAhHpKiLDAIwxHxhjbjPG/BHIFpEudRJhOJVLBKmJsQBk5RVW9o4jbMrdxOrs1ZzVpQFNPmM3pwvOfgZG3AXL37LaDbRHkVINSrASwT+BAxVsPxR4rUoiMlpE1olIhojcVcHrt4nIGhFZKSJfiEjH0MKuIf+RVUMpidYdxrv3h5YIPtn0CQ5xcEbnM2wJr8ESgZPvhvNehN0/w/MnwndTwO+PdGRKqRAESwTpxpiV5TcaY5YC6VW9UUScwLPAGKA31oilvcvt9iMw2BhzDDAd+HuIcdfMUSUCKxFkHSgI+la/8fPppk8Z2mYoLeNa2hZig3bMhXD9YugwFGbfCW+Mg9xtkY5KKRVEsEQQW8VrwSavPw7IMMZsMsYUAdOAI1pYjTFfGmMOBVa/BdKCHLN2yvUaapEQg0Ng94HgJYJlvy5j58GdnNmlkd9JXFvN2sGl78NZk2HHMnjuBPjhNW1IVqoeC5YIlojIVeU3isiVQLCZS9oB28usZwa2VeZKKumJJCJXi8hSEVmalZUV5LRVKNdryOkQWiTEhFQ19MmmT4hzxTGq/aianz9aiMCgSXDdN9C2P3x8M7x1IezfGenIlFIVCDbW0K3ADBG5hMMX/sGABzg3yHsrGomtwp+FInJp4LgjKnrdGDMFmAIwePDgmv+0LFc1BFb1ULDG4iJfEXO3zuWUDqcQ746v8emjTlJHmDgTlrwA8x6A546HMX+HYy7WIa2VqkeqLBEYY341xpwAPARsCSwPGWOGGmN2BTl2JtC+zHoacNRPQhE5FevO5bONMaG12tZUBYkgJTGG3UHaCL7a8RUHig5oI3FNOBww5Bq47mtI6QkzroFpl0De7uDvVUrViVDvI/jSGPNMYJkf4rGXAN1EpJOIeIDxwMyyO4jIAKz7E84ud+eyPcr1GgKrRBCsamj25tkkxyYzpM0QO6Nr3Fp0gctnw2mPQMbn8OwQWP1BpKNSSmHjMBHGGC9wIzAH+Bl41xjzk4g8LCIl8xg8ASQA74nIchGZWcnhwqPCqqFY9uQV4vNXXON0sPggC7Yv4PSOp0fvkBLh4nDCsJvhmv9BUjpMvxzemwQHsyMdmVJRLdT5CGrEGDMLmFVu2/1lnp9q5/mPUq7XEEBq0xj8BrIPFpbeYFbW/G3zKfQVarVQOKX2tMYq+vppWPA4bFlk9TLqqf/GSkVCwxs4rjZ83gobiwGyKulC+unmT2mX0I5+Kf1sDy+qOF1w0p1w9ZeQ0NoavG76FbB3U6QjUyrqRFkiKDqqjSAlUAqo6F6CnIIcvt35LaPTRzfe6SgjrXVfuGq+NUTF2lnwzGBrfuTc7cHfq5QKiyhMBJWUCCpoMJ63dR4+42NMpzF1El7UcnmsISpuWQHHXQUrpsEzA2HWnXAgWOc0pVRtRVci8B9dNVQ63lAFXUhnb55N52ad6Z7UvU7Ci3qJrWDM43Dzj9A/MM/B5H4w9z44uCfS0SnVaEVXIqigsTjW7aRprOuoqqFfD/7KD7/+wOhOWi1U55qlWY3HNy6BPufC4methDD/UcjPjXR0SjU60ZcInEd3AU1tGntUY/HcrXMxGEanj66r6FR5yZ3h3Ofh+m+h22nWnAf/PAYWPmHNm6yUCovoSgR+31ElAgjcVFYuEXy2+TN6JvekU7NOdRWdqkxKD7jwVbh2EaQPgy8ftRLC1/+CokNB366Uqlp0JQJfsdVtsZw2zeLIzDl8QdmRt4OVe1ZqaaC+ad0XJkyFP8yHtgNg3l/gX/2tuQ+89o5OolRjFl2JwF9cYYmga2oCv+4vZH+BNc3inC1zAPht+m/rNDwVorRB8PsPrCErWnS15j6Y3B++/Jt2O1WqBqInEfj9YPwVthF0TU0AIGN3HmBVC/Vt2Ze0RHunR1C11PEEmPQp/H6Gdbfywsfhn33hzfNhzUfgDX0uaqWima1DTNQr/sCk6o6jP3K3kkTwax4tmx/g570/c8fgO+oyOlVTItBllLXkbLXmTf7xTXh3IjRJgX4TYOBl0LJrpCNVqt6KnkTgCySCCkoE7ZPj8bgcZGTlkRvzFaDVQg1SUkc4+R4Y8WfI+AKWvWZ1Pf3mX9BxmJUQep8N7mCT6ykVXaInEZSWCI5OBE6H0LllEzb8eoCl3s8YkDqA1k1a13GAKmwcTuh+urUc+NUqJSx7HWZcbd2tfMxFMOgyq/FZKRVFbQQ+r/VYQYkAoFurRNbtzWBDzgYtDTQmia1g+G1w0zK47BMrOSx7HZ4/EaaMhKWvQMH+SEepVERFTyKooo0ArHaCbL5HEE7veHodBqbqhMMBnYbD+S/C7WutKTO9RfDJrfBUD2vWtGVvWCUIpaJM9FQNVdFGANAlpQnOxJX0SupPSnxKHQam6lx8sjV95nFXw45lsOJtWPcZrP3Eer3dIOg+2lpa99X5lVWjFz2JwB+oGqpkljF33C84Y7LoEj++DoNSESVi3ZOQNgjGPgm/rrYSwvrP4Mu/wpf/B03bQfffWkmh00na0KwapehJBKUlgoo/8srchRjjILaofx0GpeoNEevXf+u+MOJOyNsNG+bCutmw4h1rJFRXHHQeCT1GQ7ffQtM2kY5aqbCInkRQRa8hYwzzts4hztuL1dt9dRyYqpcSUmHApdbiLYQtX8H6OYESw2xrnzb9rZJC5xHWkBdaWlANVPQkgiraCFZkrWDnwZ2c0OI8vliSQ16hl4SY6PmnUUG4YqDrqdYy5u+w+2crGayfY93NvPAxqxNCm37Q/nhofxy0H6IlBtVgRM/Vroo2glmbZxHjjOHi3mOY890qvt2Yzam9W9VxgKpBEIFWva1l+O1waC9s/w62fQvbv4elL8G3z1r7Nu9gJYSSJbV3pVWTSkVS9PxVVtJGUOwr5rPNnzEibQTDurQjzr2GrzZkaSJQoYlPhh5jrAWsLqm7VlrJYft3sPkrWPWe9ZonweqR1CFQakg7FmKbRS52pQKiJxFUch/Boh2LyCnMYVzXccS4nAzpnMxXG3RaRFVDLg+kDbaWoTeAMZC7zSotbP/WSg7/e8IaABGxRk9t1dsqLaT2gtQ+kNzJujtaqToSRYmg4qqhjzd9THJsMkPbDgVgeLcUHlm3hu17D9E+Ob6uo1SNjYg1BlJSRzjmQmtb4QHY8QNs+84qPexaBWtmAsZ63RVrTcZTNjmk9oKmbfWeBmWL6EkEpUNMHP7I+wr3sWD7Ai7ucTHuQIIY0T2FR4BPV/3CtSO6RCBQ1ejFJFrdUDuPPLyt6CBkrbMaonevsZaNX8KKqYf3iW1WJjn0tpJFUicrQWgJQtVC9CSC0hLB4Y/86aZPKfYXc1aXs0q3dU1N4IQuLXj16y1cMawTHlf0jMKhIsjTBNoNtJayDu09Mjns/hlWvQ+FLx/ex+G2GqaTOkJS+tGLtkOoIKIwEVi//I0xvLPuHfq06EOv5F5H7Hr1SZ2Z9MoSZq7YyQWDdHIaFUHxydY8zenDDm8zBvbvhD3rrDkYcrZAbuBx54+Qn3PkMWKbV5AgOkJiG0hsbb2uVU5RLQoTgfWRv9/1PZv2beLRYY8i5f4TjOieQs/Wibzwv02cN6AdDof+J1H1iAg0a2ctFcnPPZwYcrYcTha7VsHaTw93nCjhirUSQmIbSGh1OEEktrFGby1Zj2mqCaORir5EEGgjmLZ2Gs1jmjO609ET1IsI143swi3TlvPcggxuHNWtLiNVqnbimltLm35Hv+b3WaWJ3G1w4Bc4sMt6zPvVev7rasj4HIryjn6vO95KCAmtoUkLiG8BccnWY+mSHFhaaOJoQKInEfgOdx/NPJDJ/O3zmdRnEjHOmAp3P7tfW+av3c0/5q1nQIckhnVtWYfBKmUThxOat7eWqhQesIbkPipZBNb3bIBD31ptGKaSYVkcrjKJIpAg4pKtJBXT1Gq7iGkKsU2PfIxJtB715rs6Ez3/0mXaCP69fDJuh5vf9fxdpbuLCH87ry9rdu7n6teX8sSF/RjbV4cMUFEiJtFags31bAwU7IND2VbbxKHswLI3sG3v4fU9GdbzglzwFQWPwd2k4iThaWIt7njrJj1P/JHPPU2s93qaBF4r2T9OSyiVsDURiMhoYDLgBF40xjxW7vUY4HVgEJANXGyM2WJLMIFEsGbfJj7d9Cl/6PsHWjWp+u7heI+LN64cwnVv/cD1by3jnP5tuW5kV3q0TrQlRKUaHJHDVVHVUVxglToK91uJpHC/NVPcUY/7Dq/n51hVWsWHrO62xYdCSyhluWKtxR1X5nmsNbKsK6bq7U6Ptbg84Iyp5HnJPjFHPne4rRKOw22Nd+Zw1aukZFsiEBEn8CxwGpAJLBGRmcaYNWV2uxLIMcZ0FZHxwOPAxbYE5Pey3yHcteSvJMcmc8Vvrgjpba2bxfLO1UP5x7z1vPbNFj5cvpOerRM5oUtLerVJpFPLJqQmxpLUxE0Tj0sblpUKhTtwoU2o5SRQvuLDSaHo4NHPy64XH4LifGs0WW++lYy8gfXifKtd5OAe8BZYS3H+4eclNQrh5CiXGJzuihNGyaPDBUNvhJ5jwx6KnSWC44AMY8wmABGZBowDyiaCccCDgefTgX+LiBhjTLiDmZGzmn+ktSUvbwdTTn+BRE/ov+o9Lgd3jenJNSd1ZvoPmXy5bjdvfreVIq//iP1EIMblINbtxO104HIIDhGcDsEhVnWTAARyhZS+L3jy0PSiVHV4AktSzQ/hovQK6TA+PBTjwovbWI+ewKObYtymGHfgucsEtlGMO/DchRen8eHCiwsfThN4xIvL+HB5vdZSdhtenPhwGh9OCnDho3BzNsf1DMM/TwUf1S7tgO1l1jOBIZXtY4zxisg+oAVwxGA/InI1cDVAhw4dahRMs2YdOSk7hQtPeZL+rY+t0TGSmni46qTOXHVSZ7w+P9v2HmLb3kPsPlDIvkPFHCgopsDrp6DYR7HP4PX58RmDMeDzG4z1OYHSwQTKPKmcCWUnpVTEeANLvs3nGd+lZte/YOxMBBX9iC1/RQtlH4wxU4ApAIMHD67RVXHUCX9i1Al/qslbK+RyOuickkDnlISwHVMppSLBzvETMoGyfdTSgJ2V7SMiLqAZsNfGmJRSSpVjZyJYAnQTkU4i4gHGAzPL7TMTuCzw/AJgvh3tA0oppSpnW9VQoM7/RmAOVvfRl40xP4nIw8BSY8xM4CXgDRHJwCoJjLcrHqWUUhWz9T4CY8wsYFa5bfeXeV4AXGhnDEoppaqmYywrpVSU00SglFJRThOBUkpFOU0ESikV5aSh9dYUkSxgaw3f3pJydy1HAf3M0UE/c3SozWfuaIypcHCnBpcIakNElhpjBkc6jrqknzk66GeODnZ9Zq0aUkqpKKeJQCmloly0JYIpkQ4gAvQzRwf9zNHBls8cVW0ESimljhZtJQKllFLlaCJQSqkoFzWJQERGi8g6EckQkbsiHY8dRKS9iHwpIj+LyE8icktge7KIzBORDYHHWszfV/+IiFNEfhSRTwLrnUTku8DnfScwDHqjISLNRWS6iKwNfNdDo+A7/mPgb3q1iEwVkdjG9j2LyMsisltEVpfZVuH3KpZ/Ba5nK0VkYG3OHRWJQEScwLPAGKA3MEFEekc2Klt4gduNMb2A44EbAp/zLuALY0w34IvAemNyC/BzmfXHgacDnzcHuDIiUdlnMvCZMaYn0A/rszfa71hE2gE3A4ONMb/BGtZ+PI3ve34VGF1uW2Xf6xigW2C5GvhPbU4cFYkAOA7IMMZsMsYUAdOAcRGOKeyMMb8YY5YFnh/AukC0w/qsrwV2ew04JzIRhp+IpAFnAC8G1gUYBUwP7NLYPm9T4CSsuTwwxhQZY3JpxN9xgAuIC8xkGA/8QiP7no0x/+PoGRor+17HAa8by7dAcxFpU9NzR0siaAdsL7OeGdjWaIlIOjAA+A5oZYz5BaxkAaRGLrKw+yfwJ8AfWG8B5BpjvIH1xvZddwaygFcC1WEvikgTGvF3bIzZATwJbMNKAPuAH2jc33OJyr7XsF7ToiURSAXbGm2/WRFJAN4HbjXG7I90PHYRkTOB3caYH8purmDXxvRdu4CBwH+MMQOAgzSiaqCKBOrFxwGdgLZAE6yqkfIa0/ccTFj/zqMlEWQC7cuspwE7IxSLrUTEjZUE3jLGfBDY/GtJsTHwuDtS8YXZMOBsEdmCVd03CquE0DxQhQCN77vOBDKNMd8F1qdjJYbG+h0DnApsNsZkGWOKgQ+AE2jc33OJyr7XsF7ToiURLAG6BXoZeLAammZGOKawC9SPvwT8bIz5R5mXZgKXBZ5fBnxU17HZwRhztzEmzRiTjvWdzjfGXAJ8CVwQ2K3RfF4AY8wuYLuI9AhsOgVYQyP9jgO2AceLSHzgb7zkMzfa77mMyr7XmcDEQO+h44F9JVVINWKMiYoFGAusBzYC90Y6Hps+44lYxcOVwPLAMhar3vwLYEPgMTnSsdrw2UcCnwSedwa+BzKA94CYSMcX5s/aH1ga+J4/BJIa+3cMPASsBVYDbwAxje17BqZitYEUY/3iv7Ky7xWraujZwPVsFVaPqhqfW4eYUEqpKBctVUNKKaUqoYlAKaWinCYCpZSKcpoIlFIqymkiUEqpKKeJQCmlopwmAhW1RKSFiCwPLLtEZEeZ9W9sON8kEckSkRer2CcucP4iEWkZ7hiUqogr+C5KNU7GmGysm7MQkQeBPGPMkzaf9h1jzI1VxJQP9A8Mm6FUndASgVIVEJG8wONIEVkoIu+KyHoReUxELhGR70VklYh0CeyXIiLvi8iSwDIshHP0CRxneWBykW52fy6lKqIlAqWC6wf0whorfhPwojHmOLFmgLsJuBVrspinjTGLRKQDMCfwnqpcC0w2xrwVGAPLadsnUKoKmgiUCm6JCQzoJSIbgbmB7auAkwPPTwV6W2OiAdBURBKNNUFQZRYD9wYm1/nAGLMh/KErFZxWDSkVXGGZ5/4y634O/5hyAEONMf0DS7sgSQBjzNvA2UA+MEdERoU5bqVCoolAoUHiuAAAAI9JREFUqfCYC5Q2AotI/2BvEJHOwCZjzL+whhU+xr7wlKqcJgKlwuNmYHCg0XcNVv1/MBcDq0VkOdATeN3OAJWqjA5DrVQdEZFJWOPGV9p9tMy+WwL77rE7LqW0RKBU3ckHxoRyQxngxmqDUMp2WiJQSqkopyUCpZSKcpoIlFIqymkiUEqpKKeJQCmlotz/A+BIKEFexkcHAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(solTchem[:,t_indx],solTchem[:,A_indx], label=\"A\")\n",
    "plt.plot(solTchem[:,t_indx],solTchem[:,B_indx], label=\"B\")\n",
    "plt.plot(solTchem[:,t_indx],solTchem[:,C_indx], label=\"C\")\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x11d7e15d0>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deZxU9Znv8c/T1dXdIN0u2IiCCBKigkSWVhONjjE6UcPV0YDBODfqGNG5GlHHzOjozRjGJBoTl0ycTHBJTC4Bo2JcrtGoATKJG40wbG5gXBplkVEEpFee+eOcaqq7q7sO0Kequ8/3/XrVq+qcOqfOcyitp3+7uTsiIpJcJcUOQEREikuJQEQk4ZQIREQSTolARCThlAhERBKutNgB7Kx9993Xhw8fXuwwRER6lUWLFn3g7tW53ut1iWD48OHU1tYWOwwRkV7FzN7u7D1VDYmIJJwSgYhIwikRiIgkXK9rIxARAWhqaqKuro76+vpih9KjVFRUMHToUNLpdORzlAhEpFeqq6ujsrKS4cOHY2bFDqdHcHc2btxIXV0dI0aMiHxebFVDZnavma03s+WdvG9m9mMzW2VmS81sQlyxAPD+Uvj+gbA2Zzgi0svU19czcOBAJYEsZsbAgQN3upQUZxvBL4BTunj/VGBU+JgG/DTGWGDuRdDwMTx0YayXEZHCURLoaFf+TWKrGnL3P5rZ8C4OOQP4pQfzYL9gZnuZ2f7u/n63BnLDnm23N7y6Y98Nm7r1UiIivVExew0NAd7N2q4L93VgZtPMrNbMajds2LBzV7n4P2HPA8msuuAAew2DS/688xGLiGRZu3YtU6dOZeTIkYwePZrTTjuN119/HYDbbruNiooKNm3a8Qfn/PnzMTPuueee1n2LFy/GzPjhD38IwPnnn8+IESMYN24cEyZM4Pnnn4/9PoqZCHKVX3KukuPuM929xt1rqqtzjpDu3P6fgXT/8HPCfen+MPjwnfscEZEs7s6ZZ57JCSecwOrVq1m5ciXf+973WLduHQCzZ8/myCOP5OGHH25z3tixY7n//vtbt+fMmcMRRxzR5phbbrmFJUuWcNNNN3HxxRfHfi/FTAR1wIFZ20OB92K5Uv0mPhowkkubLqd54CGw7aNYLiMiyTFv3jzS6TSXXHJJ675x48Zx3HHHsXr1arZs2cKNN97I7Nmz25w3bNgw6uvrWbduHe7Ok08+yamnnprzGscffzyrVq2K9T6guN1HHwUuM7M5wNHApm5vH8i4+jWeqX2XJx5cyrXnXsuB+/SP5TIiUhzfeWwFK9/7uFs/c/QBVfzL/xrT6fvLly9n4sSJOd+bPXs255xzDscddxyvvfYa69evZ9CgQa3vT548mQceeIDx48czYcIEysvLc37OY489xtixY3fvRiKIs/vobOB54BAzqzOzC83sEjPLpM8ngDeBVcBdwP+JKxaAstLgVhtbtsd5GRER5syZw9SpUykpKeGss87igQceaPP+2WefzQMPPNCaMNr71re+xbhx45g5c2ab9oS4xNlrqOPdtX3fgUvjun57ZakwETQrEYj0NV395R6XMWPG8OCDD3bYv3TpUt544w1OPvlkABobGzn44IO59NIdP3eDBw8mnU7z9NNPc8cdd/Dcc8+1+YxbbrmFyZMnx3sDWRIz11BriUCJQES6wYknnkhDQwN33XVX676FCxcyffp0brjhBt566y3eeust3nvvPdasWcPbb7edBXrGjBncfPPNpFKpQofeQWKmmEiHJYImVQ2JSDcwMx5++GGuuOIKbrrpJioqKhg+fDjz58/npz9tOz72zDPPZM6cORx99NGt+4455phCh9wpc8/ZY7PHqqmp8V1ZmGbR8ldo/s35lEz5BUeOPSyGyESkkF555RUOO0z/L+eS69/GzBa5e02u4xNTNTRs+U840l7jgP/6cbFDERHpUfp+1dCNg6C5gWoAgyGrfg03/BpKy+H69cWOTkSk6Pp+iWD6Ujh8CttLKwBoTlXA2CkwfVmRAxMR6Rn6fiKoHAzllVhzA/WeJtXSAOVVULlfsSMTEekR+n4iANi6nq2fOY8zG2ew6qCzYcu6YkckItJj9P02AoCps9i2uYFXXnqGFw49jVGfG17siEREeoxklAjYMaCsQQPKRKQbbNy4kXHjxjFu3DgGDx7MkCFDWrffeeednNNTz58/n0mTJhU79A6SUSJgxxQTTS29a9yEiHSjzWvhwQtg8i92u51w4MCBLFmyBIAbbriBAQMGcPXVV+PuHHPMMZx33nnMmTMHgCVLlrROT90TJa5EoCkmRBJswQ/gnRdgwc2xXaKr6akBtmzZwuTJkzn00EM599xzyQzqXbRoEX/1V3/FxIkT+dKXvsT77weTMZ9wwglceeWVHH/88Rx22GEsXLiQs846i1GjRnH99dd3S8yJKRGkSoxUiWmKCZEkCscTtaq9J3jEMJ6oq+mpIViRbMWKFRxwwAEce+yx/PnPf+boo4/mm9/8Jo888gjV1dXcf//9XHfdddx7770AlJWV8cc//pE77riDM844g0WLFrHPPvswcuRIrrzySgYOHLhbMScmEQCkU6ZpqEWSaPpSeOp6ePVxaN4Gpf3gsEnw198teChHHXUUQ4cOBYKSwltvvcVee+3F8uXLW2csbWlpYf/992895/TTTweC1c3GjBnT+t7BBx/Mu+++q0SwM8pSJaoaEkmicDwRLQ1QWhE8xzSeqLPpqTOyF6FJpVI0Nzfj7owZM6bT9Ykz55SUlLQ5v6SkhObm5t2OOTFtBBC0E6hEIJJQW9fDxAvgG88EzzGNJ+pseuoFCxZ0es4hhxzChg0bWhNBU1MTK1asiCW+XFQiEJFkmDprx+tJt8Z2mc6mp7799ttZs2ZNznPKysp48MEHufzyy9m0aRPNzc1cccUVjBlTmAV3EjMNNcAJt8zjM0P34sfnjO/mqESk0DQNdec0DXUX0qkS9RoSEWknUYmgrFRVQyIi7SUvEahEINJn9Laq7ULYlX+TRCWCaj7k2rVXweaeO9RbRKKpqKhg48aNSgZZ3J2NGzdSUVGxU+clqtfQ2Vtnc2jTimB4eYy9BkQkfkOHDqWuro4NGzYUO5QepaKionXAWlTJ6DXUfnh5hparFJGE6KrXUJclAjObEOHzm9y9Z6/7GA4vb1jxKOXeUNTh5SIiPU2+qqEFwELAujhmBDC8uwKKRTi8PO2NNJCmXMtVioi0ypcIFrr7iV0dYGZ/6MZ44rN1Pc/vfQY/23IcvzxipZarFBEJdZkI8iWBqMf0CFNn8fjcZaxcuQ4mfaPY0YiI9Bh5u4+a2WAzGxy+rjazs8ysMBNgdLPy0hIam1uKHYaISI/SZSIws4uB54EXzOzvgceBScBcM7uwAPF1Kw0oExHpKF8bwWXAGKAf8DbwKXdfa2Z7A/OAe2KOr1ulU6Y1i0VE2slXNdTk7p+4+0ZgtbuvBXD3D4G8v6hmdoqZvWZmq8zsmhzvDzOzeWa22MyWmtlpu3QXEZWlUrRsd1q2KxmIiGTkSwTbzSwdvv5yZqeZVeQ718xSwJ3AqcBo4BwzG93usOuB37j7eGAq8O87EftO0wL2IiId5UsEZxH+5e/udVn7BwL/kOfco4BV7v6muzcCc4Az2h3jQFX4ek/gvShB76p0KhgOoXYCEZEd8nUffSd728yqwnO2AS/n+ewhwLtZ23XA0e2OuQH4vZl9E9gDOCnXB5nZNGAawLBhw/JctnPlKhGIiHQQafZRM7vYzNYBS4FF4SPfhD+5RiO3r5w/B/iFuw8FTgN+ZWYdYnL3me5e4+411dXVUULOqbVqSCUCEZFWUWcfvRoY4+4f7MRn1wEHZm0PpWPVz4XAKQDu/nzY9rAvEMtMcOlUkAiaVCIQEWkVdT2C1cAnO/nZC4FRZjbCzMoIGoMfbXfMO8AXAczsMKACiG1OWZUIREQ6iloiuBZ4zsxeBFrnc3b3yzs7wd2bzewy4CkgBdzr7ivMbAZQ6+6PEjQ432VmVxJUG53vMc6LXZZSG4GISHtRE8HPgD8Ay4DIv6Lu/gTwRLt93856vRI4Nurn7a60SgQiIh1ETQTN7n5VrJEUQLlKBCIiHURtI5hnZtPMbH8z2yfziDWyGGRKBE0qEYiItIpaIvha+Hxt1j4HDu7ecOKlNgIRkY7yLVW5v7u/7+4jChVQnDTFhIhIR/lKBPeGM43OB54E/uTuzbFHFZPMOAI1FouI7NBlG4G7nwqcQJAIziRYl2Bu2F6w63M9FEm/hg3cXzYD2xLLeDURkV4pb2Oxu9e7+5PuPt3dawj6/pcCPzGzl2KPsBvtvfB2jrTXOPS1WCc5FRHpVWx3xm+ZWVk4s2jB1NTUeG1tvmmO2rlxEDQ3dNxfWg7Xq3QgIn2fmS0K/5jvIN+aApvN7OPwsTlre7OZfVzoJLDLpi+Fw6fgpf0AaCqpgLFTYPqyIgcmIlJ8+aahrixUILGqHAzlldDSQL2nKd/eAOVVULlfsSMTESm6qOMIMLMjgOPCzT+6+9J4QorJ1vUw8QLO/POn+P5Bixi3ZV2xIxIR6REiJQIzmw5cBMwNd80ys5nu/m+xRdbdps7CgNXP/44nDzqRcaceWuyIRER6hKglgguBo919K4CZ3Qw8D/SeRBAqKy3RgDIRkSxR5xoyoCVru4XcK5D1eOmUaa4hEZEsUUsEPwdeNLOHw+2/Ae6JJ6R4qUQgItJWpETg7rea2Xzg8wQlgQvcfXGcgcWlrLREU0yIiGSJ3GsI+AvQHJ5jZjbB3V+OJ6z4pFNKBCIi2aL2GvpX4HyCtYszQ5EdODGesOJTllLVkIhItqglgrOBkb1mJHEXytVGICLSRtReQ8uBveIMpFDSqRL1GhIRyRK1RPB9YLGZLQdaZ29z99NjiSpG6jUkItJW1ERwH3AzsAzo1b+iZaUlbGnotWvriIh0u6iJ4AN3/3GskRRIWo3FIiJtRE0Ei8zs+8CjtK0a6nXdRzWOQESkraiJYHz4/Nmsfb2y+2i5SgQiIm1EHVn8hbgDKRT1GhIRaSvfCmWT8n1AlGN6EvUaEhFpK1+J4BYzW0PXM41+D3i8+0KKlxKBiEhb+RLBOuDWPMe80U2xFERQNeT5DxQRSYh8axafUKA4CibTa8jdMeuVSyqIiHSrqFNM7BIzO8XMXjOzVWZ2TSfHnG1mK81shZn9Os54AMpSwY+/upCKiAR2ZhrqnWJmKeBO4GSgDlhoZo+6+8qsY0YB1wLHuvuHZjYorngyykqD3NfU4pTHdvciIr3HLpcIzCyd55CjgFXu/mY4a+kc4Ix2x1wE3OnuHwK4+/pdjSeqslRwy2owFhEJ7FQisMCJZnY3wV/5XRkCvJu1XRfuy/Zp4NNm9mcze8HMTunkutPMrNbMajds2LAzIXeQbi0RKBGIiEDERGBmR5vZHcDbBNNM/CdwaL7Tcuxr312nFBgFnACcA9xtZh2mu3b3me5e4+411dXVUULulEoEIiJt5RtQ9l0ze4NgrMAygqkmNrj7fZnqnC7UAQdmbQ8F3stxzCPu3uTufwFeI0gMscm0ETQoEYiIAPlLBNMIxhL8FPh/7r6Rjn/Vd2YhMMrMRphZGTCVoDSR7bfAFwDMbF+CqqI3I37+LsmUCFQ1JCISyJcIBgPfBU4HVpnZr4B+Zpa3v427NwOXAU8BrwC/cfcVZjbDzDIL2jwFbDSzlcA84FthsolNpkSgqiERkUC+AWUtwO+A35lZBTAJ6A+sMbNn3f1rec5/Anii3b5vZ7124KrwURB7NH3A/WUz2L75l/SR1TdFRHZL5F5D7l7v7g+6+1cI6vGfii+s+AxffidH2msMXnxHsUMREekRLPijvPeoqanx2tranT/xxkHQ3NBxf2k5XB/78AURkaIys0XuXpPrvVinmOhRpi+Fw6ewPVUBQEuqAsZOgenLihyYiEhx5es+OiV8HlGYcGJUORjKK7GWBuo9TUlLA5RXQeV+xY5MRKSo8pUIrg2fH4o7kILYup7Nh3+dMxtn8Obwr8KWdcWOSESk6PJ1A91oZvOAEWbWfgwA7n56jnN6rqmz2PLRNl6p/QMLR09i5FHDih2RiEjR5UsEXwYmAL8CfhR/OPFLa0CZiEgb+cYRNAIvmNkx7r7BzCqD3b6lMOF1P00xISLSVtReQ/uZ2WJgObDSzBaZ2eExxhWb8szIYpUIRESA6IlgJnCVux/k7sOAfwj39TqtVUPNvWv8hIhIXKImgj3cfV5mw93nA3vEElHMUiVGqsRobGkpdigiIj1C1MUa3zSz/0vQaAzwt8Bf4gkpfmWpEk06JyISiloi+DugGpgbPvYFLogrqLilU0ZTi6qGREQgYokgXITm8phjKZiy0pR6DYmIhJIz11CW8lJVDYmIZCQyEQRVQ0oEIiKQ0ERQphKBiEirSG0EZlYNXAQMzz7H3f8unrDilU6VaECZiEgoavfRR4D/BJ4Ben0H/LLSElUNiYiEoiaC/u7+T7FGUkBlqRL1GhIRCUVtI3jczE6LNZICUhuBiMgOURPBdIJkUG9mm8PHx3EGFqeylKqGREQyog4oq4w7kEJSiUBEZIeobQSY2enA8eHmfHd/PJ6Q4qdeQyIiO0SqGjKzmwiqh1aGj+nhvl6prLSEJpUIRESA6CWC04Bx7r4dwMzuAxYD18QVWJzKSlUiEBHJ2JmRxXtlvd6zuwMpJE1DLSKyQ9QSwfeBxWY2DzCCtoJrY4sqZioRiIjsELXX0Gwzmw8cSZAI/snd18YZWJxUIhAR2aHLqiEzOzR8ngDsD9QB7wIHhPt6pXSqhO0OLdu1OI2ISL4SwVXANOBHOd5z4MRuj6gAykqD/NfYvJ1+ZakiRyMiUlxdJgJ3nxa+PNXd67PfM7OKfB9uZqcAdwAp4G53z9nl1MwmAw8AR7p7bZTAd4cSgYjIDlF7DT0XcV8rM0sBdwKnAqOBc8xsdI7jKgmWwXwxYiy7rTxMBNuaev1EqiIiuy1fG8FgM5sI9DOz8WY2IXycAPTP89lHAavc/U13bwTmAGfkOO5fgR8A9Tnei0VVvzQAm+ubCnVJEZEeK18bwZeA84GhwK1Z+zcD/5zn3CEEDcsZdcDR2QeY2XjgQHd/3Myu7uyDzGwaQVsFw4YNy3PZ/PYME8HHSgQiInnbCO4D7jOzr7j7Qzv52ZbrI1vfNCsBbiNINF1y95nATICamprd7upTVVFKNR8y8vGz4euzoHK/3f1IEZFeK+o4gofM7MvAGKAia/+MLk6rAw7M2h4KvJe1XQkcDsw3M4DBwKNmdnrcDcZV/dJcXvowe26ohQU3w6Rb858kItJHRV2z+D8I2gS+ANwNTAZeynPaQmCUmY0A1gBTga9l3nT3TcC+WdeYD1wde6+hGwcxsrmBkZk7r70neJSWw/XrY720iEhPFLXX0DHu/nXgQ3f/DvA52v6134G7NwOXAU8BrwC/cfcVZjYjnNK6OKYvpWXMZLZ5WbBd2g/GToHpy4oWkohIMUWdayjTo+cTMzsA2AiMyHeSuz8BPNFu37c7OfaEiLHsnsrBpCqqKKeJJisj3dIA5VVqJxCRxIqaCB4zs72AW4CXCRp974otqrhtXc/c1F9TN2IqV+z9J9iyrtgRiYgUTd5EEPbuedbdPwIeMrPHgYqwjr93mjqLn926gE/ZAJh0VrGjEREpqrxtBOFiND/K2m7o1UkgVNUvzaZtGkcgIhK1sfj3ZvYVC/t59gV79ktrQJmICNHbCK4C9gCazayeYLCYu3tVbJHFrKqilFXrm4sdhohI0UUdUFYZdyCFpqohEZFApKohM3s2yr7eZM9+aTbXN7Fdi9OISMJ1WSII1xzoD+xrZnuzY/6gKuCAmGOLVVVFmu0OWxubqaxIFzscEZGiyVc1dDFwBcGP/iJ2JIKPCdYa6LWq+gW3vmlbkxKBiCRavtlH7wDuMLNvuvu/FSimgmidinpbM+xd5GBERIooamPxv5nZMcDw7HPc/ZcxxRW7qgqtSSAiAtFnH/0VMBJYAmTWd3Sg9yaC1hKBEoGIJFvUcQQ1wGh37zNdbDJVQ+pCKiJJF3Vk8XKChWP6jB1VQxpUJiLJFrVEsC+w0sxeAhoyO929eOsK7KYBFcGtq2pIRJIuaiK4Ic4giiFVYlRWlKpqSEQSL2qvoQVmdhAwyt2fMbP+QCre0OJXVaGJ50REok4xcRHwIPCzcNcQ4LdxBVUoVf3SwTgCEZEEi9pYfClwLMGIYtz9DWBQXEEVyp79StVGICKJFzURNLh7Y2bDzEoJxhH0aqoaEhGJnggWmNk/A/3M7GTgAeCx+MIqjKBqSIlARJItaiK4BtgALCOYiO4J4Pq4giqUPbUmgYhI5O6j/YB73f0uADNLhfs+iSuwQqiqSLO1sYXmlu2UpqLmRBGRviXqr9+zBD/8Gf2AZ7o/nMLKTEW9WaOLRSTBoiaCCnffktkIX/ePJ6TC0XxDIiLRE8FWM5uQ2TCzicC2eEIqnH39Q+4vm8En/72m2KGIiBRN1DaCK4AHzOy9cHt/4KvxhFQ4h7z+H1Tba6xdeDt8+t+LHY6ISFFY1JmlzSwNHEKwXOWr7l6U+pSamhqvra3dvQ+5cRA0N3TcX1oO16/fvc8WEemBzGyRu9fkem9nusocCXwGGA+cY2Zf747gimL6Ujh8CttLg/bv5pIKGDsFpi8rcmAiIoWXzBXKKgdDeSXW0kC9pynf3gDlVVC5X7EjExEpuFhXKDOzU4A7CGYqvdvdb2r3/lXAN4BmggFrf+fub+/MNXbZ1vUw8QImPz+Kfx26kPFb1hXksiIiPU3URJBZoez9qB8cDjq7EzgZqAMWmtmj7r4y67DFQI27f2Jmfw/8gEI1Qk+dhQHvLX6ah/b/POP/ZmxBLisi0tPEuULZUcAqd38TwMzmAGcArYnA3edlHf8C8LcR4+k2VRWlmopaRBItzhXKhgDvZm3XAUd3cfyFwO9yvWFm04BpAMOGDduFUDq3Zz/NQCoiyRap15C7LwBeBSrDxyvhvq5Yro/KeaDZ3xK0Q9zSyfVnunuNu9dUV1dHCTmyKk08JyIJF3WFsrOBl4ApwNnAi2Y2Oc9pdcCBWdtDgffaH2RmJwHXAae7e47O/fGqqtBU1CKSbFGrhq4DjnT39QBmVk0w6dyDXZyzEBhlZiOANcBU4GvZB5jZeILlL0/JfHahVfVL87EmnRORBIs6oKyk3Q/1xnznunszcBnwFPAK8Bt3X2FmM8ws08h8CzCAYPqKJWb26M6Fv/uq+pWqakhEEi1qieBJM3sKmB1uf5VOGnazufsTBIvYZO/7dtbrkyJePzZVFWkam7dT39RCRTpV7HBERAouUiJw92+Z2VnA5wkagWe6+8OxRlYgmamoP97WpEQgIonUZfWOmX3KzI4FcPe57n6Vu18JbDSzkQWJMGaDKssBeH9TfZEjEREpjnxtBLcDm3Ps/yR8r9cbOWgAAKs3bMlzpIhI35QvEQx396Xtd7p7LTA8logKbNg+/SktMSUCEUmsfImgoov3+nXxXq+RTpUwfu96zlh8EWzWxHMikjz5EsFCM7uo/U4zuxBYFE9IhXdpyVw+Vb8MFtxc7FBERAquyxXKzGw/4GGgkR0//DVAGXCmu6+NPcJ2umWFsgytVCYiCbHLK5S5+zp3Pwb4DvBW+PiOu3+uGEmg24UrlTWnghqw7aX9tFKZiCRO1HEE84B5eQ/sbcKVylKZlcqa67VSmYgkzs6sWdw3bV1P4/gLOLNxBisOmAxaqUxEEibqFBN919RZlAMblz/DfXsfyy1Tjih2RCIiBaUSQWhk9QCNJRCRRFIiCI0ctAerN2ylq15UIiJ9kRJBaGT1ADZta2Lj1sZihyIiUlBKBKGR1eGcQ+tVPSQiyaJEEBo5aADVfMjwx6ZoqgkRSRQlgtD+VRVcVfZbqj98WVNNiEiiqPsowI2DKGlu4JxMWqy9J3hoqgkRSQCVCKB1qolGCxapQVNNiEiCKBFA61QTaW+k3tN4S4OmmhCRxFAiyNi6nrcPnsqZjTNY/+mvaaoJEUkMJYKMqbPY96s/4e30wdxtX4Ft/63eQyKSCEoEWQaUl/Llsftz8Cv/jr/zgnoPiUgiqNdQthsHcUtzAxjgqPeQiCSCSgTZpi/FD59MPeo9JCLJoUSQrXIwVl5FOY3UeynevA1KStV7SET6NCWC9rauZ9sR5/HM9onB9tvPFTceEZGYqY2gvVVP07+5gUmpcPujt+GGPdVOICJ9lkoE7YWjjFvCBe2bKMP77wvf+EORAxMRiYcSQXvhKOOS7Y00WRml3giffADP/QR+fqrGFohIn6OqoVy2rscw0t4YdCUFWDo7eP7RITBkIpjBl2+DJ/8RJv9CDcoi0mvFWiIws1PM7DUzW2Vm1+R4v9zM7g/ff9HMhscZT2RTZ8FVK+HwKXRcuNJhTS3ULYS534B3XoCn/yUoLby/DO76Itx9UvA6sy/X8+Z1sHlt21JG++2MzvbnsjPHdvf5vfHcYl67N55bzGsn/Z539x66EFsiMLMUcCdwKjAaOMfMRrc77ELgQ3f/FHAb0HOG8oZVRIbhVpIjIQAbXgXfHpQW3n4O/9nnOyaJzp4X3AwLfrDjNXTczuhsfy47c2x3n98bzy3mtXvjucW8dtLveXfvoQsW12LtZvY54AZ3/1K4fS2Au38/65inwmOeN7NSYC1Q7V0EVVNT47W1tbHE3MGcc2HAfkEbwcpHWpNBprbIPaghiotnXStbA2kmVT3UZt/jH3+FcpoiHZvL7pzfG88t5rV747nFvLbuuRM72ZPRzBa5e03O92JMBJOBU9z9G+H2/waOdvfLso5ZHh5TF26vDo/5oN1nTQOmAQwbNmzi22+/HUvMnWqXEIAOJYTMrBStSaLdvvbPLRhgpNhOA2VsLt2bqpYPKfNGGqyc/9rjOJ7ZeyonfTibI7b+iXJvCPd/nrnVf8/Hpfu0uX5V80bO2vAfkY7NZXfO743n9ta49e+VzHtuCStvUmwPZjw4bBL89Xd3qm2yq0QQZ2Nxrj9mc/1+5jsGd58JzISgRLD7oe2kqbOC5znnwl4HwZAJ2DsvQMNmKKuELe8DbW/Gus26KuYAAAbCSURBVHq2FClvARxKKyhvaaR8r71h44bW7aMOPYijJk2Bx56Dl+dl7R/OUZNOzh3nY7+Lfmx3n98bz+2tcevfK3H3nGquD/aVVkAM66XEmQjqgAOztocC73VyTF1YNbQn8N8xxrR7MgkhW6a08MEb8N7LQWIwgiSRKocB1bBlQ9vnPQYFA9Uyn1n7c3j1/8PEC6DmgmA7sx7C1vW59+eyM8d29/m98dzeGrf+vZJ3z3PODfZlfi+6eb2UOKuGSoHXgS8Ca4CFwNfcfUXWMZcCY939EjObCpzl7md39bkFbSMQEekjilI15O7NZnYZ8BSQAu519xVmNgOodfdHgXuAX5nZKoKSwNS44hERkdxiHVDm7k8AT7Tb9+2s1/XAlDhjEBGRrmmKCRGRhFMiEBFJOCUCEZGEUyIQEUm42LqPxsXMNgC7OrR4X+CDvEf1LbrnZNA9J8Pu3PNB7l6d641elwh2h5nVdtaPtq/SPSeD7jkZ4rpnVQ2JiCScEoGISMIlLRHMLHYARaB7TgbdczLEcs+JaiMQEZGOklYiEBGRdpQIREQSLjGJwMxOMbPXzGyVmV1T7HjiYGYHmtk8M3vFzFaY2fRw/z5m9rSZvRE+713sWLuTmaXMbLGZPR5ujzCzF8P7vd/MyoodY3cys73M7EEzezX8rj+XgO/4yvC/6eVmNtvMKvra92xm95rZ+nDlxsy+nN+rBX4c/p4tNbMJu3PtRCQCM0sBdwKnAqOBc8xsdHGjikUz8A/ufhjwWeDS8D6vAZ5191HAs+F2XzIdeCVr+2bgtvB+PwQuLEpU8bkDeNLdDwWOILj3Pvsdm9kQ4HKgxt0PJ5jWfip973v+BXBKu32dfa+nAqPCxzTgp7tz4UQkAuAoYJW7v+nujcAc4Iwix9Tt3P19d385fL2Z4AdiCMG93hcedh/wN8WJsPuZ2VDgy8Dd4bYBJwIPhof0tfutAo4nWMsDd29094/ow99xqBToFy541R94nz72Pbv7H+m4QmNn3+sZwC898AKwl5ntv6vXTkoiGAK8m7VdF+7rs8xsODAeeBHYz93fhyBZAIOKF1m3ux34R2B7uD0Q+Mjdm8PtvvZdHwxsAH4eVofdbWZ70Ie/Y3dfA/wQeIcgAWwCFtG3v+eMzr7Xbv1NS0oisBz7+my/WTMbADwEXOHuHxc7nriY2SRgvbsvyt6d49C+9F2XAhOAn7r7eGArfagaKJewXvwMYARwALAHQdVIe33pe86nW/87T0oiqAMOzNoeCrxXpFhiZWZpgiQwy93nhrvXZYqN4fP6YsXXzY4FTjeztwiq+04kKCHsFVYhQN/7ruuAOnd/Mdx+kCAx9NXvGOAk4C/uvsHdm4C5wDH07e85o7PvtVt/05KSCBYCo8JeBmUEDU2PFjmmbhfWj98DvOLut2a99ShwXvj6POCRQscWB3e/1t2Huvtwgu/0D+5+LjAPmBwe1mfuF8Dd1wLvmtkh4a4vAivpo99x6B3gs2bWP/xvPHPPffZ7ztLZ9/oo8PWw99BngU2ZKqRd4u6JeACnAa8Dq4Hrih1PTPf4eYLi4VJgSfg4jaDe/FngjfB5n2LHGsO9nwA8Hr4+GHgJWAU8AJQXO75uvtdxQG34Pf8W2Luvf8fAd4BXgeXAr4DyvvY9A7MJ2kCaCP7iv7Cz75WgaujO8PdsGUGPql2+tqaYEBFJuKRUDYmISCeUCEREEk6JQEQk4ZQIREQSTolARCThlAhERBJOiUASy8wGmtmS8LHWzNZkbT8Xw/XON7MNZnZ3F8f0C6/faGb7dncMIrmU5j9EpG9y940Eg7MwsxuALe7+w5gve7+7X9ZFTNuAceG0GSIFoRKBSA5mtiV8PsHMFpjZb8zsdTO7yczONbOXzGyZmY0Mj6s2s4fMbGH4ODbCNcaEn7MkXFxkVNz3JZKLSgQi+R0BHEYwV/ybwN3ufpQFK8B9E7iCYLGY29z9T2Y2DHgqPKcrlwB3uPuscA6sVGx3INIFJQKR/BZ6OKGXma0Gfh/uXwZ8IXx9EjA6mBMNgCozq/RggaDOPA9cFy6uM9fd3+j+0EXyU9WQSH4NWa+3Z21vZ8cfUyXA59x9XPgYkicJ4O6/Bk4HtgFPmdmJ3Ry3SCRKBCLd4/dAayOwmY3Ld4KZHQy86e4/JphW+DPxhSfSOSUCke5xOVATNvquJKj/z+erwHIzWwIcCvwyzgBFOqNpqEUKxMzOJ5g3vtPuo1nHvhUe+0HccYmoRCBSONuAU6MMKAPSBG0QIrFTiUBEJOFUIhARSTglAhGRhFMiEBFJOCUCEZGE+x9e1VE9kDYjdwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,1],label='CAMP')\n",
    "n=1000\n",
    "plt.plot(solTchem[:n,t_indx],solTchem[:n,A_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of A [mol/m3]')\n",
    "plt.legend(loc='best')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x11d9ca3d0>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXhU9dXA8e+ZJZMACWELq6yi7GvccK07ijtaLFZEq/ZtVdzaaqu+VqnVulVbu1B3XwruimilbrgvBMEIKArIvgRQlgQy63n/mJsQQpKZQG4mkzmf57lP5t65d3ImA/fMbxdVxRhjTObypDoAY4wxqWWJwBhjMpwlAmOMyXCWCIwxJsNZIjDGmAznS3UA9dW+fXvt2bNnqsMwxpi0Mnfu3E2q2qGm59IuEfTs2ZOioqJUh2GMMWlFRFbU9pxVDRljTIazRGCMMRnOEoExxmS4tGsjMMYYgHA4zOrVqykvL091KE1KdnY23bp1w+/3J32NJQJjTFpavXo1ubm59OzZExFJdThNgqqyefNmVq9eTa9evZK+zqqGarJ9PTw2GrZvSHUkxphalJeX065dO0sCVYgI7dq1q3cpyRJBTd79E6z8BN69K9WRGGPqYElgT3vzN7GqoapuL4BocNd+0SPxzReAm0pSF5cxxrjISgRVfN/rVBSI4t11sG1vmPRlymIyxjRd69evZ9y4cfTp04cBAwZwyimn8M033wBw//33k52dzdatWyvPnz17NiLCI488Unls3rx5iAj33HMPABdddBG9evVi2LBhjBgxgo8//tj192GJAGByAdzamrZLXkAAL9Fdz32/DO49IH6OMcY4VJWzzjqLY445hqVLl7Jo0SLuuOMONmyIty1OmzaNgw46iBdffHG36wYPHszTTz9duT99+nSGDh262zl333038+fP58477+Tyyy93/b1YIgCYVAyDziVIAAAFYs6fJoZAi/bws7dTGKAxpql555138Pv9/PznP688NmzYMI488kiWLl1KaWkpkydPZtq0abtd1717d8rLy9mwYQOqyuuvv87o0aNr/B1HHXUUS5YscfV9gLURxOV2IpLVCr+GCHuy8GsIIUZI4o91xyak6FEYc1+qIzXG1OD3ryxk0dptDfqaA7rk8b+nDaz1+QULFjBy5Mgan5s2bRrnn38+Rx55JIsXL6akpISCgl21CmPHjuXZZ59l+PDhjBgxgkAgUOPrvPLKKwwePHjf3kgSrETg2L5pLf8XPY7PjnsW8nsAQpaGEEAg3mh8a2urIjLGJDR9+nTGjRuHx+Ph7LPP5tlnn93t+fPOO49nn322MmFU96tf/Yphw4YxZcqU3doT3OJqiUBETgYeALzAw6p6Z7XnuwNPAPnOOTeo6mtuxlSbp/v8kTu/+ZqiYYfBEcXxsQSzbiK4YAYBgqgvB+k/Bk78QyrCM8bUoa5v7m4ZOHAgzz333B7Hi4uL+fbbbznhhBMACIVC9O7dm1/+8peV53Tq1Am/388bb7zBAw88wEcffbTba9x9992MHTvW3TdQhWslAhHxAg8Bo4EBwPkiMqDaaTcBz6jqcGAc8De34qnT9vX86KMJjGwXon0rp4iW2wkCuWQRolz9ECmHQB7kdkxJiMaYpuXYY48lGAzyr3/9q/LYnDlzmDRpErfeeivLly9n+fLlrF27ljVr1rBixe6zQN92223cddddeL3e6i/d6NysGjoYWKKqy1Q1BEwHzqh2jgJ5zuPWwFoX46mVzr6LvsEFXB94afcnykqg8GKuzr2HN/w/Qhe+aKONjTFAfODWiy++yBtvvEGfPn0YOHAgt956K7Nnz+ass87a7dyzzjqL6dOn73Zs1KhRnHnmmY0Zcq1EVd15YZGxwMmq+jNn/6fAIap6RZVzOgP/BdoALYHjVXVuDa91GXAZQPfu3UdWz6x7bXIBRIJ7Hq82gOzZolWUv3Q1F/jeQgovtkZjY5qAr776iv79+6c6jCappr+NiMxV1cKaznezRFDTOOfqWed84HFV7QacAjwlInvEpKpTVLVQVQs7dKhxpbW943QbjXiyAYj5cmDwubsPIJtcwLkzB/FT35sIao3Gxphmx81EsBrYr8p+N/as+rkEeAZAVT8GsoH2Lsa0O6cdwBMLEsSPRIN7tgM4ySIkzhiDmpKFMcakMTcTwRygr4j0EpEs4o3BM6qdsxI4DkBE+hNPBBtdjGlPZSW86D2Je3v8DRk5EUqrtQE4ycKvFY3GO+G79xs1RGOMcZNriUBVI8AVwCzgK+K9gxaKyG0icrpz2nXApSLyBTANuEjdarSoxZqTH+a6sgvpcsBB8br/cVP3PMlpNP6Z/07W+HpC6XqbmdQY02y4Oo7AGRPwWrVjt1R5vAg43M0YEpm38gcACnu2rf2kcVORyQX8X8RmJjXGND8ZP7J40/b4zb1Lfk7dJ04qZuv+Z7JTs+L71lZgjGkmMj4RbCuPAJCbnaBwlNuJ1q3bEpAwIbLi6xbYADNjMtbmzZsZNmwYw4YNo1OnTnTt2rVyf+XKlTVOTz179mzGjBmT6tD3kPGTzm3bGaZFlhe/N4mcWFbCoi5jueO7/XmyzaP4tjTQeAZjTOPYvh6emwhjH9/nL3Ht2rVj/vz5ANx66620atWK66+/HlVl1KhRTJgwoXIQ2fz58yunp26KrERQHqZ1jj+5k8dNpe15f2G0dw7eHSXO5HTGmLTRCMvQ1jU9NUBpaSljx46lX79+jB8/nor+MXPnzuXoo49m5MiRnHTSSaxbtw6AY445hmuuuYajjjqK/v37M2fOHM4++2z69u3LTTfd1CAxW4lgZ4S87CQTweQCukSC/LTir2YNxsakh+qzCLj4f7eu6akhviLZwoUL6dKlC4cffjgffvghhxxyCFdeeSUvv/wyHTp04Omnn+Z3v/sdjz76KABZWVm89957PPDAA5xxxhnMnTuXtm3b0qdPH6655hratWu3TzFnfCLYujNMXk6Sf4ZJxTDrJiKLXsEXKyfmy8bT/zSbkdSYps75v8vXM+NjgXw5kKLZhA8++GC6desGxEsKy5cvJz8/nwULFlTOWBqNRuncuXPlNaefHu9xP3jwYAYOHFj5XO/evVm1apUlgn21rTxMp7zs5E52Bpd5Y/HBZVkRazA2Ji04/3eJBsGX7Wpnj9qmp65QdREar9dLJBJBVRk4cGCt6xNXXOPxeHa73uPxEIlE9jlmayMoD5OXbBsBQFkJUjiRWwoeYIb/ZLT6SGRjTNNUVgIjJ8LP3oz/dOn/bm3TU7/77ru1XnPggQeycePGykQQDodZuHChK/HVJLMTwfb13L/jt3Txbk3+mnFTYcx9DBl5BH/YPoYdWzfa1NTGpAPn/y6dBtc+i0ADqG166i5dutR6TVZWFs899xy/+c1vGDp0KMOGDdtjsRo3uTYNtVsKCwu1qKioQV5LX7kGLXqM4k5nM+x/Hq3XtZtLg/znrp/wE+9beGxqamManU1DXbv6TkOdmW0ETg8CAURg2Ibn4dbnk+9BMLmAdpEgF1QsLGS9h4wxaSwzq4acqaVjvvi0EhFPdv2mi6i+joG3ntcbY0wTkpmJwOlBIJFyytWPN1bPHgTVeg/ZdBPGpEa6VW03hr35m2RmIgAoK2HDgeM5K3Qb6w/4Sf17EDi9h/7Y9a/8R4629YyNaWTZ2dls3rzZkkEVqsrmzZvJzk6yS7wjM9sIAMZNpXjher76Yi6bjx5P566t6309wPB5a/j++WdBf4gPW7dGY2MaRbdu3Vi9ejUbNzbuWlZNXXZ2duWAtWRlbiJg18yjSU8xUd3kAs6MBHf9Fa3R2JhG4/f76dWrV6rDaBbqTAQiMiKJ1wiralq2km7bGQZIftK56pxh68GFMwhoEPXlICkatm6MMXsrUYngXeJrD0sd5/QCejZUQI1pq5MIWiVai6A2TqNxlrOecSBSbo3Gxpi0k+gOOEdVj63rBBF5uwHjaVTbysPkBnx4PXXluQTKSoiOvJifzDmAG9u8w0ELX4Sjf2PJwBiTNursNZQoCSR7TlO1bWekfvMM1WTcVHyn3UePgYeybEsMLd9iC9sbY9JKwjoREekEoKrrRaQDcCSwWFUbb0Ykl2wrDydeojIZkwu4v2Kuc8UajY0xaaXOEoGIXA58DHwiIv8DzATGAC+IyCWNEJ+rtu2s58yjtZlUTHTgWHZiC9sbY9JPoq/DVwADgRxgBbC/UzJoA7wDPOJyfK7aVh6ha37Ovr9Qbie82XkECMcbjaNBxBqNjTFpItHI4rCq7lDVzcBSVV0PoKo/EK8ESWvbdtZjveJEykrYdOB4Lg5dT3lWO7CF7Y0xaSJRiSAmIn5VDQOnVhwUkWyawfQU2+qzTGUi46bSQZWxS8YRCG6yhe2NMWkj0V3wbJxv/qq6usrxdsB1bgXVGKIxZXuwHgvXJzK5AIkEObti3xqMjTFpIlH30ZWqWrkgpojkiUhbYCfwudvBuam0YnqJhqoaqpzaOj7ZU72ntjbGmBRJql7E6T10G/EEUNE2oEBvl+Jy3bby+KjivIboPgqVo4w90RBB/PjrO7W1McakSLL1/NcDA1W1p6r2cra0TQKwa3qJBisRQOXi2K8e/BTPR44gsuAFm5raGNPkJZsIlgI73Aykse0qETRgInAWxz78iB9RTgCPjTI2xqSBZOtFbgQ+EpFPgWDFQVW9ypWoGsE+zzxam8kFdIwE+alNTW2MSRPJlgj+CbwNfALMrbKlrW07KxqLG3hJhor1jL3xRuOorWdsjGnikr0LRlT1WlcjaWSVVUMNXSKoup4xfrKi5fDd+w37O4wxpgElWyJ4R0QuE5HOItK2YnM1Mpdt2xlGBFplubBIW1kJMnIiUw6YwhLthpaut7YCY0yTlexd8CfOzxurHEvz7qMRcgM+PPuyFkFtxk2FyQVcFXlk15I+1lZgjGmiEs0+2hmgSpfRXs2l+2iDzTxam0nF6KCxlBOI79uMpMaYJipRieBRZ6bR2cDrwAdVRxqns23l4YbtOlpdbickkEeAXctY2oykxpimKNEUE6OBY4gngrOIr0vwgtNe0N398NyztSFnHq1NWQk7hkzg7PBtLGw/Gha+aAPMjDFNTsLGYlUtV9XXVXWSqhYSn2zOB/xVRD5zPUKXxJepdKGhuKpxU2l59gMU9C1k8fdRW8bSGNMk1XsqaVX9TlX/pqqnA0fUda6InCwii0VkiYjcUMs554nIIhFZKCL/rm88e8u/cwO/Xnet+9/QJxfw+IoTOCc2C9FYvMH41tYwucDd32uMMUlK1Fi8XUS2Odv2KvvbRWSbqobquNYLPASMBgYA54vIgGrn9CXeE+lwVR0IXL3P7yhJ44NP02tHsfvf0CcVExtky1gaY5quOutGVDV3H177YGCJqi4DEJHpwBnAoirnXAo85Kx4hqq6369ycgFEgpzfWN06czvhCVRZxjJiy1gaY5qWpKuGRGSoiFzhbEOSuKQrsKrK/mrnWFUHAAeIyIci8omInFzL775MRIpEpGjjxo3JhlyzScUE+5/DTm3Eb+hlJZQOupCzQrfxdcdTrNHYGNOkJJUIRGQSMBUocLapInJlostqOFZ9nWMf0Jd4z6TzgYdFJH+Pi1SnqGqhqhZ26NAhmZBrl9uJoLcFAcJEPQGINsK6AeOmkjf2QfJ7jeDrTRFrNDbGNCnJdpu5BDhEVcsAROQu4GPgL3VcsxrYr8p+N2BtDed84qyJ/J2ILCaeGOYkGddeiW0v4f+ix3HASVdy6PczoLQRvp1PLmBaJLhr30YaG2OaiGSrhgSIVtmPUvM3/qrmAH1FpJeIZAHjgBnVznkJ+BGAiLQnXlW0LMmY9tqCI/7GLZGL8XQeAmPui08J4bZJxUQGWqOxMabpSbZE8BjwqYi86OyfCTxS1wWqGhGRK4BZgBd4VFUXishtQJGqznCeO1FEFhFPLr9S1c1780bqY9fMoy6PI6gqtxO+7Dw81mhsjGlikroTqup9IjKb+LgBASaq6rwkrnsNeK3asVuqPFbgWmdrNBWL0uS6OcVETcpKKB18IePm9uOO7nMZ1hhVUsYYk0B9vhJ/B0Sca0RERqjq5+6E5a6d4XgtVwu/t3F/8bip5AEF2z7jsrUH8OGlx9LIqcgYY/aQbK+h24Fi4EHgXme7x8W4XBWOxgDw++o9sLpB/PTQHpRsD/LGIisRGGNSL9kSwXlAn7pGEqeTcDTei9XvdWEtgiQcc2ABXfNzmPHBXE4pegjGPm5tBcaYlEn2K/ECYI/+/emqskTgSU2JwOsRxh/ancPXPIau+NjGFBhjUirZEsEfgXkisgCo7AzvTDyXdsLRGF6PuLM6WTImF/CLSHDXX9/GFBhjUijZRPAEcBfwJRBzL5zGEY5qyqqFAJhUDLNuIrRwBlkaRH05SP8xcOIfUheTMSZjJZsINqnqg65G0ojC0Rh+b2qqhQDI7QSBXPy6a/Uy16e5MMaYWiSbCOaKyB+JjwyuWjWUlt1Hw9EYWalMBABlJUjhxdyyYiSHfD+Ds7ZvqP/iEMYY0wCSTQTDnZ+HVjmmwLENG07jCEcUXyqrhqByWovjF67nsqfyyRk0glNSG5ExJkMlO7L4R24H0phSXjVUxXH9O9K9bQse/eA7ThncOdXhGGMyUKIVysYkeoFkzmlqwjFNfdWQw+sRLhrVk6IVP7Bw8WJ4bLStVWCMaVSJ7oZ3i8hwERlR2wbc0RiBNqRwJJb6qqEqzi3sRquAj+9fmwwrP7FxBcaYRpWoamgDcF+Cc75toFgaTVOqGgLIvacrCyQIW50DNq7AGNOIEq1ZfEwjxdGoQk0sETCpmJ0zb4SvZ5IjofhaBTauwBjTSJrQ3bDxRKJNp40AgNxO5LTKJyDxtQo00gjLZxpjjKMJ3Q0bTzjatNoIACgrYdvAn3JW6Da+6HR24yyfaYwx1G89gmYjHI3RMtDE3vq4qeQDXXcUMXHF/nx48bG0SHVMxpiMkLBEICI9nPWEEZFDReR6ETnL/dDcE59rqGkWhn5+dG9+2BHmmTmrUh2KMSZDJBpHcDPwNvCJiEwG/gy0B64SkT83QnyuCEdjZPmaWNWQo7BnWwp7tGHKe8sIRdJ+fj9jTBpI9LX4fKA/UAhcBRyrqjcAJwDHuxyba8LRGL4UrUWQjCuO3Z+1W8t54fPVqQ7FGJMBEt0Ny1U1pKpbgKWqugNAVSNA2q5W1pSrhgCOPqADQ7q15m+zlxKJWqnAGOOuRHfDfBE5W0TOAfKcxxX7rRshPlc05aohABHhymP7svL7Hbw8f22qwzHGNHOJus68C5zmPH6vyuOK/bTU1KuGAI7vX0D/znk89M4SzhzeFW+qVlMzxjR7iUYWT2ysQBpTU68agopSwf78YurnvDnnC05adKMtcm+McUXTvhu6JBSN4W/CVUMVTh7Yib4FrQi+eSdqk9EZY1zSxEZVNY5IU1ihLAmeOzryRiS464BNRmeMcUGicQTnOj97NU447ovGlJjS5NsIAJhUjA4aSzkBANSXA4PPhUlfpjgwY0xzkuhueKPz83m3A2ksYac7ZjpUDZHbCQnkESC+yD22yL0xxgWJqoY2i8g7QC8RmVH9SVU93Z2w3FORCNKhagiAshIovJgblg7nyG0zOXP7BrypjskY06wkSgSnAiOAp4B73Q/HfeGoAjT5XkOVxk1FgLHfbuKCR9pS2mMgE1IdkzGmWUnUfTREfJ6hUaq6UURy44e1tHHCa3gVJYImNw11Aofv345DerXlL28vYezIbk1v9lRjTNpK9mtxRxGZBywAFonIXBEZ5GJcrolsXcfTWbeRG/4+1aHUi4hww+h+bCoN8q/3l6U6HGNMM5JsIpgCXKuqPVS1O3Cdcyzt5H56HwfJYgYt+UeqQ6m34d3bcMrgTkx5bxkl28tTHY4xpplINhG0VNV3KnZUdTbQ0pWI3DK5AG5tTd6CJ/GI0nv5dLi1dfx4GvnVSf0IRWI8+Na3qQ7FGNNMJJsIlonIzSLS09luAr5zM7AGN6kYBp1LzJsNQNSbnZZ98nu1b8n5B3dn2merWLoxbZtqjDFNSLKJ4GKgA/CCs7UH0mseotxOEMhFokHK1Y8nmr4LxF91XF+yfR7u+s/XqQ7FGNMMJNX1RFV/IL4wTXorK6HkwPFcVDyQRwYtoEuaLhDfITfAz4/uw71vfMNHSzcxqk/7VIdkjEljadKZvoGMm8rSg37PV9qDlYfdDuOmpjqivXbpUb3pmp/D7TO/IhrTVIdjjEljriYCETlZRBaLyBIRuaGO88aKiIpIoZvxQHzmUUijAWW1yPZ7ufGUfny1bhtP20L3xph94NrdUES8wEPAaGAAcL6IDKjhvFzi1U6fuhVLVRFnZHHaTDFRh1MHd+bgnm2597+L2VYeTnU4xpg0ldTdUEQ6iMhvRWSKiDxasSW47GBgiaouc0YoTwfOqOG824E/AY3SMT6tJp1LQES45bQBfL8jxINvWndSY8zeSfZr8cvE1yh+E3i1ylaXrkDVOovVzrFKIjIc2E9VZ9b1QiJymYgUiUjRxo0bkwy5ZhVVQ2kxDXUSBnVtzY8L9+Pxj5bzzYbtqQ7HGJOGkr0btlDV36jqM6r6fMWW4JqavnJXtmqKiAe4n/go5Tqp6hRVLVTVwg4dOiQZcs3CzahqqMKvT+5Hq2wfN7+0AFVrODbG1E+yd8OZInJKPV97NbBflf1uwNoq+7nAIGC2iCwHDgVmuN1gHGlGVUMV2rbM4tcn9ePT777npflrUh2OMSbNJJsIJhFPBuUist3ZtiW4Zg7QV0R6iUgWMA6oXNNAVbeqantV7amqPYFPgNNVtWgv3kfSws2saqjCuIP2Y+h++fxj5sdEHjkZtqfnGAljTONL6m6oqrmq6lHVbOdxrqrmJbgmAlwBzAK+Ap5R1YUicpuIpGxBm1AzrBoC8HiEyWcM4oLgdDyrbKF7Y0zykp7U3rl5H+Xszk7UwAugqq8Br1U7dkst5x6TbCz7ojlWDQEwuYDBkSCDKz5RW+jeGJOkZLuP3km8emiRs01yjqWdcDMZULYHZ1I99eUAUE6A6KD0m1TPGNP4ki0RnAIMU9UYgIg8AcwDah0t3FRVVA35PM2sRFBlUr2oJ0BWNMSCzcqQNJxUzxjTuOrztTi/yuPWDR1IYwlHY/i9gkgzSwQQX+h+5ES8l73FB/mns37NCpaU2FTVxpi6JVsi+CMwT0TeIT4+4CjgRteiclEkGmt+1UIVqkyi1/9n/+L4+96l7/PFPHP5YXiaWwnIGNNgku01NI14P/+K9QgOU9XpbgbmlnBUm28iqKJDboBbxgygaMUPPPpheq0hZIxpXHXeEUWkn/NzBNCZ+CCxVUAX51jaCTlVQ5ng7BFdOb5/AXfPWmyrmRljapXoq/G1zs97a9jucTEu14QjzbhqqBoR4Y6zB5OT5eW6Z76o7DprjDFV1XlHVNXLnIejVfVHVTfiPYnSTiSWGVVDFQpys7ntjEHMX7WFKe8vS3U4xpgmKNk74kdJHmvyQtEYvgypGqpw2pDOnDq4M/e/8Q0L1mxNdTjGmCYmURtBJxEZCeSIyHARGeFsxwAtGiXCBhaOxJrd9BKJiAiTzxxE25ZZTJo+j52haKpDMsY0IYnuiCcRbwvoBtzHrvaBa4HfuhuaOzKtaqhCm5ZZ3HfeMJZtKuP2VxelOhxjTBNS5zgCVX0CeEJEzkli/YG0EM6gXkPVHb5/ey47sjf/fG8ZRx/QgZMGdkp1SMaYJiCpAWWq+ryInAoMBLKrHL/NrcDcEorE8GVgiaDCdSceyIdLN/Gb54sZ3LU1XfJzUh2SMSbFkp107h/Aj4EriY8sPhfo4WJcrglHM6+NoKosn4cHxw0nHIlxxb8/r5yEzxiTuZK9I45S1QuBH1T198Bh7L76WNqItxFkZtVQhd4dWvHHc4bw+cot3D1rcarDMcakWLKJoNz5uUNEugBhoJc7IbkrlEEDyupy+tAujD+kO1PeW8abi2w1M2MyWbJ3xFdEJB+4G/gcWA5McysoN4Wb86Rz9XTzmAEM7JLHtc/MZ+XmHakOxxiTIgnviCLiAd5S1S1Oz6EeQL/aVhpr6uKTzmV21VCFbL+Xv48fCcBlTxWxIxRJcUTGmFRImAicxWjurbIfVNW0HZ7arKeh3gvd27XgwfOHs3jDdn7z/JfotnXw2GjYbtVFxmSKZO+I/xWRc6QZrOYSimpGdx+tyTEHFnD9iQfyyhdr+Wr6zbDyE3j3rlSHZYxpJMkuTHMt0BKIiEg58S6kqqp5rkXmknj30bTPZw3uFx+O4pfZQVjrHCh6JL75AnBTSUpjM8a4K9mFaXJV1aOqWaqa5+ynXRIAqxqqjUwqJjzgHMoJABDzZcPgc2HSlymOzBjjtmQHlL2VzLEmb/t6HuN/yY/9kOpImp7cTvhzWhMgRBA/RIKUe1pCbsdUR2aMcVmi2UezRaQt0F5E2ohIW2frCXRpjAAbkr57F4WymCPWPJLqUJqmshKk8GK+O3MG02LHM//rbygP20ylxjR3oqq1PykyCbia+E1/DfG2AYBtwL9U9a+uR1hNYWGhFhUV1e+iyQUQCe553Oq/azWzeC1X/Hsepw7pzF/GDcfjsXYVY9KZiMxV1cKanku0QtkDqtoLuF5Ve6tqL2cbmooksNcmFcOgc1FffIK1sMfqvxMZM6QLN47ux6vF6/jDa1+lOhxjjIuSnX30LyIyCuhZ9RpVfdKluBpWbicI5MbrvdVPIBaEQJ7Vfydw2VG9Wbe1nEc++I7OrbP52ZG9Ux2SMcYFSSUCEXkK6APMByoqjRVIj0QAUFZC+dAJnP3ZAdy//zz6ldqAqUREhJvHDGDDtnImv/oV7VsFOHN411SHZYxpYMmOIygEBmhdDQpN3bipbNm6k68+fZt5g8fQ7+DuqY4oLXg9wv0/HsaWHXO47tkvaBnwccIAK0kZ05wk26F+AZD2y1lFovE8ZuMI6ifb7+VfEwoZ1LU1v/z353y0ZFOqQzLGNKBk74jtgUUiMktEZlRsbgbmhpCzCItNOld/rQI+nph4EL3ateRnTxYxd8X3qQ7JGNNAkq0autXNIBpLuDIRWIlgb+S3yHzTF9kAABfPSURBVOKpSw7mx1M+YcKjc3jykoMZ0b1NqsMyxuyjZKeYeJf4GgR+5/Ec4usSpBWrGtp3BXnZTLv0UNq1ymLCI58xf9WWVIdkjNlHyU4xcSnwHPBP51BX4CW3gnKLVQ01jE6t48mgTcssfvrIp8xbaVN2GJPOkv1q/EvgcOIjilHVb4ECt4JySzhiVUMNpUt+DtMuO5S2LbO44OFP+WTZZti+3tYyMCYNJXtHDKpqqGJHRHzExxGklbBVDTWorvk5PHP5YXTOz+Gixz5jzcu/t7UMjElDyTYWvysivwVyROQE4BfAK+6F5Y5wzKqGGlrHvGzeKDsX8QZhiXPQ1jIwJq0k+9X4BmAj8CVwOfAacJNbQbnFqobcIVcXE+q/ay2DiM3lZExaSfaOmAM8qqrnqupY4FHnWJ1E5GQRWSwiS0Tkhhqev1ZEFolIsYi8JSI96hd+/VjVkEtyO5HVojUBCRMiC080yPyNMbRV2jUjGZORkr0jvsXuN/4c4M26LhARL/AQMBoYAJwvIgOqnTYPKFTVIcR7Jf0pyXj2Sth6DbmnrAQZORHPpW/ycdszWL9mJb9+rpiQUwozxjRdybYRZKtqacWOqpaKSIsE1xwMLFHVZQAiMh04A1hU5XXeqXL+J8AFScazV2xAmYvGTQXi/6BGXfU497/5Lc++9S1rtuzk7+NH0rqFP7XxGWNqlewdsUxERlTsiMhIYGeCa7oCq6rsr3aO1eYS4D9JxrNXrGqocYgI155wAPeeO5Q5y7/n7L9/yIrNZakOyxhTi2TviFcDz4rI+yLyPvA0cEWCa2qqf6mxy6mIXEB8htO7a3n+MhEpEpGijRs3JhnynqxqqHGdM7IbT11yCJtKQ5zx0Id8aJPVGdMkJTvFxBygH/A/xLuO9lfVuQkuWw3sV2W/G7C2+kkicjzwO+B0Va1hPUlQ1SmqWqiqhR06dEgm5BpVJgKflQgay6G92zHjisMpyA1w4aOf8cgH35HOs5kb0xzV5454EDAEGE684ffCBOfPAfqKSC8RyQLGAbvNWCoiw4lPW3G6qrre4byiaijLqoYaVY92LXnhF4dzXL8Cbp+5iGuens+OUCTVYRljHK6tUKaqERG5ApgFeIl3P10oIrcBRao6g3hVUCvi1U4AK1X19L19M4lUlAh8thB7o2sV8PGPC0by13eWcP+b37Bo3Tb+fsFI+nRolerQjMl4rq5QpqqvER98VvXYLVUeH1+f19tX4WgMkfiqW6bxeTzCVcf1ZXj3fCZNn8/pf/mAO84ezBnDusbnKXpuIox93NaSNqaRZdQKZeGo4vd6cEofJkWO7NuBmVceQb/OeUyaPp9fPfsF4XfutHmKjEmRZEsEFSuUfQZUNui6WY3jhnA0Zu0DTUSX/ByevuxQmNwR38LQridsniJjGl3GrVDms66jTYbP64Frv2Tj87+i1fLXySFE2BPAO+B0PCf9IdXhGZMx6rNC2ddArrN95RxLKxVVQ6YJye1Eh3btySZCSLLwRkPMWlLGsvKWqY7MmIyR7Apl5wGfAecC5wGfishYNwNzQ3Z5CVMiN9vCKU1NWQlSOJGsy99mRe9x+Ms3MfqB95ny3lKiMRtzYIzbJJmOQCLyBXBCRV9/EekAvKmqQ12Obw+FhYVaVFS0V9d+8OefMmrLK3gKL4Yx9zVwZKahbNhWzk0vLeCNRRsY2q01d40dQr9OeakOy5i0JiJzVbWwxueSTARfqurgKvse4IuqxxrLXiWCyQUQqWHQsjVINlmqyszidfzvjIVs2xnmkiN7Mem4vrTISrZZyxhTVV2JINkK89dFZJaIXCQiFwGv4vIEcQ1qUjEMOpegxBdOwZdjC6c0cSLCaUO78Na1R3P2iK78891lnHDfe7yxaINNUWFMA0u2sfhXxKeCGAIMBaao6q/dDKxB5XaCQC5+DRHCD9EgBPJs4FIaaNMyiz+NHcozlx9Giywvlz5ZxITH5rCkpDTxxcaYpNSZCERkfxE5HEBVX1DVa1X1GmCziPRplAgbSlkJb7U8lRva/hlGToRSazBOJwf3astrk47k5jEDmLfiB07+83vcPnMRW3aE4qOSHxttnQCM2Ut1thGIyEzgt6paXO14IfC/qnqay/HtYV8ai8/7x8d4PDD9ssMaOCrTmDaVBrn79cU8M3cVedl+pnV5hv5rn0dGTrROAMbUoq42gkQtbz2rJwEAVS0SkZ4NEFujCsditPJbY2O6a98qwF1jh/DHr47HEw3CGucJG5VszF5J1EaQXcdzCRevb2psionmxXN1vBNA1Bv/Z7pTs3jTdzRvn/wmMRt/YEzSEt0V54jIpdUPisglQKKFaZqccERtionmxOkE4I2FUF+AbAlTJjlc/NwqTnnwfV5fsM4SgjFJSFRPcjXwooiMZ9eNvxDIAs5yMzA3hGMxm2KiuSkrgZETkcKJUPQYp23fQKzfUP7y1hJ+/n+fc2DHXP7nmD6MGdI5PreRMWYPyQ4o+xEwyNldqKpvuxpVHfalsfjIP73NQT3act+PhzVwVKapicaUGV+s4W/vLOXbklL2a5vDZUf25pyR3WxQmslI+zygTFXfUdW/OFvKksC+sqqhzOH1CGcN78asq4/iXxcW0r5VgJtfXsioO9/m7llfs2Fb+a6TrfupyXAZ9dUoHLWqoUzj8QgnDOjI8f0LKFrxAw+/v4y/zV7KP99dxujBnZlwWA9GLrgLqVgUx7qfmgxkicBkBBHhoJ5tOahnW5ZvKuPJj1fwm6KjCCwO7zrJup+aDJVRd8VwVMnyZdRbNjXo2b4lt5w2gOhV8/mu8ymUE5+DaqdmMS//BL44532bz8hklIy6K4ajMXy2cL1xtGjXjV5dOhGQMDFvgICE+fp7OOOJJRx/37v8ffZS1m3dmeowjXFdxiQCVSUSsxXKTDVlJcjIiXgufQtP4cWM7ZfFn84ZQtuWWdz1+teMuvNtxk35mOmfrYzPa1SdNTSbZiBj2gjC0XhR36qGzG7GTd31eMx9+IkvwXfeQfuxYnMZL81by0vz13DDC19y88sLOGL/9owZ0oXjB3SkdY4f3v0TWEOzSXMZlAhiAPit+6hJUo92LZl0fF+uOm5/vlyzlZnF63i1eB3XPfsFXwcmgFhDs2keMi4R+DxWIjD1IyIM6ZbPkG753HByP+av3sKUz1+k/5d3cUTkU3IkRDkBVhQcCydO5gBVROwLh0kfGZQI4lVDfqsaMvvA4xFGdG/DiO5Hot7nYe4HRCSLLA3x2dowNz+8mI55y/nRgQUcdUAHRvVpR36LrFSHbUydMiYRRLet4+ms21gf/jvQI9XhmGZAyjZC4cX4nHmOzt2yjkC/IcxeXMKrxeuYPmcVHoHB3fI5Yv92jOrTnpE92pDt98ZfYPt6eG4ijH3cVsszKZXUXENNyd7ONbTtuStp9eVTfNfzPPpMnOJCZMbsEo7G+GLVFt77dhPvf7uR4tVbicbi41iG75fPIb3acv7GB+i0ZJotqGMaRV1zDTX/RDC5ACLBPY9bo55pRKXBCJ99t5mPlmzm10VHkUV4j3OingAbrlpB59bZ1sZgGlxmJ4Lt62HWTcS+egVPtJyoNxvvgNPgxD9Ycdykxvb1hP7zW7yLX8UbLSdIgFmxQm4PjWcj+XTMCzBsv3jj9NBu+Qzu2prWLfx7vIZVK5n62JelKtOfs3iJRIOUq59ANAiBPPvPY1IntxNZOa0hFgJfNoFoiFMPOpAew05l3sof+HzlFopXb2HWwl2D1Lq3bcHgrq0Z0CWPAV3yOGThH8lZ+Qli4xdMA2j+iQCgrISNB45nQvFAHh60gK6lNgrUpJizoA5OQ7O3dAND98tn6H75XHR4/JQtO0IUr97Kl2u2snDtVorXbOHeb04ku4bxC1FPFvMnfEPfjq3Iy/bX/DsrWGnCVNP8q4YcHy/dzPn/+oR/X3oIo/q0dyEyY9y3beMqQq/9lvwV/8UXK6ecAP+NHcTtoZ+wkXwAOuYF2L+gFft3aEWfglb0bt+KXh1a0jkvG49HYOa1MPexeCKy0kTGyOyqIUckFh9QZovXm3SW12E/aNselserlbKjIU476ACGjTqLbzZs55uS7SwpKWVpSSnPf76G0mCk8tqvAxNqLE3EvAEiN65PfvoVK1E0OxmTCCpHFlsiMOmuWrWSlG6ge7sWdG/XguMH7Loxqyol24Ms21jGsk2l/HPtCxyy5H6GlX1INkF2ahavRw/ijvLxbLr5P3TKy6Zbmxy65ufQtU0OXfNb0Dk/my6tc+iSn01uRZXTvsyvZEmkScqYRBCKOCOLba4hk+6qTZRXGxGhY142HfOyOaxPO6AHvPIifD4b9QbIjoQZNbAXvzngaFZ9v4NVP+xg9Q87mbP8B14pXkc0tnu18eLABAI1tk8EmD/hawpys2nfKkBOlrf22G2SviYpYxKBVQ0ZQ2VpQpzSRMfSDYwd2W2P0yLRGCXbg6zbupM1W8pZt2Unf938PEcue4AhpR/sXqKIjGfj3z+uvLZVwEf7Vlm0bxWgXass2rUKcFvxsfi0yjTeThJRXwBJdjzPvpQmrCRSp4xJBLtmH7VEYDJYkqUJn9dDl/wcuuTnMLJyRpY+8MpL8Pk74I23T5w8rC8HHDKGku1BNm4LsrE0yKbSIBu3B9lcGuK7TWXMWf4Db5bfz2+9UznJW0SOhHarltp+839o0yKL/BZZ5Of4yW/hp3VOfMur8nPkl7fTZcXHbH39dspPvIfcbB8tsrzJDb7b15JIM09CmZMInKohn1UNGbP3qrRPSNFj5JRuYGCX1gxMcFk0poRf/ozAFx8T8wbIjoYY2KsrF/U+hC07QvywI8yWHWG27gyxpKSUrTvDbN0ZJhiJ7dHInb/wKVj4FOXqp0/oCVoFfORm+2kV8NEq20fLgI+WWV5aBnz8cdHx+GsoiUQ9AT4b/xUtA15aZHnJ9nvJ8XvJyfKS7fPGe1dVtS+JpKHaVFDXEoqr3UdF5GTgAcALPKyqd1Z7PgA8CYwENgM/VtXldb3m3nYffXXWLI78aALhC1+jXZ8R9b7eGLOPpo+HVh0rG7kp3bB7CaUG5eEopZtWE3j7Flosex1vtJyIN5tVBcfyXq+r2aj5bC8PUxqMUhoMUxqMUBqMUhaMUBaMkBPcyKToE5zoqVYSiYyv7G5bk2y/h2y/l09jPyFQw3QgIcni+gP+S7bfQ8DnJeDzEPB7yPJ6nZ8eLnz7YHyxPVe1i3kDzLtwMVleD1k+D36v4Hce+zyC3+fB74kf9/7nOmTu4/HkC/vU7TclU0yIiBf4BjgBWA3MAc5X1UVVzvkFMERVfy4i44CzVPXHdb3u3iaCH+4eTn7pMqLtD8R35Wf1vt4Yk0KvXAOfPw7eLIiG6nUz1Feuhs+fAI8foiG2DLiA5Yfezs5QlJ3hKDucnztD8cfl4fi2MxzFv2MDJ699iCGlHxDQIEECfJx1GP/Mvph10TzKwzFC0RjBcJTySGy3BvYO/MDvfDVUhyVIQlBDV9+a1HO+tFSNIzgYWKKqy5wgpgNnAIuqnHMGcKvz+DngryIi2pDZ6dbWALQBEPBtXlx5jFu3NtivMca4qFqXWeoxO4CUbdzt2jalG2jTvU2SVw+CV2bG20Wc6UCOGdKHY8acUePZ0ZgSisQIRWIEo1FyZn1K9oJd1WGjBvbigYNPqjwnElPC0RjBSIxINP44HI3x9I6ZHLLkfvpsfgd/LEjMWV7eQwx8OdB/THy+tAbiZiLoCqyqsr8aOKS2c1Q1IiJbgXbApqonichlwGUA3bt3r18Ul78P03+Cbl2FAApIfncYN61+r2OMSZ0kG7kb/FqoVxLyeoScLK/ThdYPke+h8OLdeml1TGpmg94Q6gKbwuDLxhMpjx/2ZYML86W5mQhqapWt/k0/mXNQ1SnAFIhXDdUris5DwN+i8hcJgL8FdBpUr5cxxmSoVCWhqglo+vhdr1fPElEy3EwEq4H9qux3A9bWcs5qEfEBrYHvGzyS8q3QoR8c/Zt4y/3OLQ3+K4wxpkFVTSJXF+967MJAPDcTwRygr4j0AtYA44CfVDtnBjAB+BgYC7zdoO0DFa5fvOvxoLMb/OWNMSaduZYInDr/K4BZxLuPPqqqC0XkNqBIVWcAjwBPicgS4iWBcW7FY4wxpmauDihT1deA16odu6XK43LgXDdjMMYYUzebb8EYYzKcJQJjjMlwlgiMMSbDWSIwxpgMl3ZrFovIRmDFXl7enmqjljOAvefMYO85M+zLe+6hqh1qeiLtEsG+EJGi2iZdaq7sPWcGe8+Zwa33bFVDxhiT4SwRGGNMhsu0RDAl1QGkgL3nzGDvOTO48p4zqo3AGGPMnjKtRGCMMaYaSwTGGJPhMiYRiMjJIrJYRJaIyA2pjscNIrKfiLwjIl+JyEIRmeQcbysib4jIt87PZNfpSwsi4hWReSIy09nvJSKfOu/3aRHJSnWMDUlE8kXkORH52vmsD8uAz/ga59/0AhGZJiLZze1zFpFHRaRERBZUOVbj5ypxDzr3s2IRGbEvvzsjEoGIeIGHgNHAAOB8ERmQ2qhcEQGuU9X+wKHAL533eQPwlqr2Bd5y9puTScBXVfbvAu533u8PwCUpico9DwCvq2o/YCjx995sP2MR6QpcBRSq6iDi09qPo/l9zo8DJ1c7VtvnOhro62yXAX/fl1+cEYkAOBhYoqrLVDUETAdqXn06janqOlX93Hm8nfgNoivx9/qEc9oTwJmpibDhiUg34FTgYWdfgGOB55xTmtv7zQOOIr6WB6oaUtUtNOPP2OEDcpyVDFsA62hmn7OqvseeKzTW9rmeATypcZ8A+SLSeW9/d6Ykgq7Aqir7q51jzZaI9ASGA58CHVV1HcSTBVCQusga3J+BXwMxZ78dsEVVI85+c/usewMbgcec6rCHRaQlzfgzVtU1wD3ASuIJYCswl+b9OVeo7XNt0HtapiQCqeFYs+03KyKtgOeBq1V1W6rjcYuIjAFKVHVu1cM1nNqcPmsfMAL4u6oOB8poRtVANXHqxc8AegFdgJbEq0aqa06fcyIN+u88UxLBamC/KvvdgLUpisVVIuInngSmquoLzuENFcVG52dJquJrYIcDp4vIcuLVfccSLyHkO1UI0Pw+69XAalX91Nl/jnhiaK6fMcDxwHequlFVw8ALwCia9+dcobbPtUHvaZmSCOYAfZ1eBlnEG5pmpDimBufUjz8CfKWq91V5agYwwXk8AXi5sWNzg6reqKrdVLUn8c/0bVUdD7wDjHVOazbvF0BV1wOrRORA59BxwCKa6WfsWAkcKiItnH/jFe+52X7OVdT2uc4ALnR6Dx0KbK2oQtorqpoRG3AK8A2wFPhdquNx6T0eQbx4WAzMd7ZTiNebvwV86/xsm+pYXXjvxwAznce9gc+AJcCzQCDV8TXwex0GFDmf80tAm+b+GQO/B74GFgBPAYHm9jkD04i3gYSJf+O/pLbPlXjV0EPO/exL4j2q9vp32xQTxhiT4TKlasgYY0wtLBEYY0yGs0RgjDEZzhKBMcZkOEsExhiT4SwRGGNMhrNEYDKWiLQTkfnOtl5E1lTZ/8iF33eRiGwUkYfrOCfH+f0hEWnf0DEYUxNf4lOMaZ5UdTPxwVmIyK1Aqare4/KvfVpVr6gjpp3AMGfaDGMahZUIjKmBiJQ6P48RkXdF5BkR+UZE7hSR8SLymYh8KSJ9nPM6iMjzIjLH2Q5P4ncMdF5nvrO4SF+335cxNbESgTGJDQX6E58rfhnwsKoeLPEV4K4Eria+WMz9qvqBiHQHZjnX1OXnwAOqOtWZA8vr2jswpg6WCIxJbI46E3qJyFLgv87xL4EfOY+PBwbE50QDIE9EcjW+QFBtPgZ+5yyu84KqftvwoRuTmFUNGZNYsMrjWJX9GLu+THmAw1R1mLN1TZAEUNV/A6cDO4FZInJsA8dtTFIsERjTMP4LVDYCi8iwRBeISG9gmao+SHxa4SHuhWdM7SwRGNMwrgIKnUbfRcTr/xP5MbBAROYD/YAn3QzQmNrYNNTGNBIRuYj4vPG1dh+tcu5y59xNbsdljJUIjGk8O4HRyQwoA/zE2yCMcZ2VCIwxJsNZicAYYzKcJQJjjMlwlgiMMSbDWSIwxpgM9/8mVa3tImTWmwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,3],label='CAMP')\n",
    "plt.plot(solTchem[:n,t_indx],solTchem[:n,B_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of B [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x11daabe50>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXhU5dn48e89M8lMgLAmIUCAgIAKAoIBFJcqFgVLpVrcWvsWN/T3VgWttVrbV6rUDatSX19bV3ADK26IVhEEbBWRBDDssigQEAh7IOtM7t8fcwghJJkBMplk5v5c11wzZ5tzn5xc557zPM95HlFVjDHGxC9XtAMwxhgTXZYIjDEmzlkiMMaYOGeJwBhj4pwlAmOMiXOeaAdwrFJSUjQzMzPaYRhjTKOSk5OzU1VTq1vW6BJBZmYm2dnZ0Q7DGGMaFRHZWNMyKxoyxpg4Z4nAGGPinCUCY4yJc42ujqA6ZWVl5OXlUVxcHO1QGhSfz0dGRgYJCQnRDsUY04BFLBGIyEvACGCHqp5WzXIBJgGXAIXAaFVdfDz7ysvLIzk5mczMTIJfa1SVXbt2kZeXR5cuXaIdjjGmAYtk0dBkYFgty4cD3Z3XGODZ491RcXExbdq0sSRQiYjQpk0bu0sypjEr2AYvD4eC7Ud+rmMRSwSq+jmwu5ZVRgKvaNBXQEsRaXe8+7MkcDT7mxhTR07kInwi285/DDZ9BfMfPfJzHYtmHUEHYHOl6Txn3g9VVxSRMQTvGujUqVO9BGeMaWAKtsH062DUZEhuW3/bwpEX4RFPHLVYVQmUK2UBpTRQTqm/HH95OWV+pcVnf6H5xgXs+vABNp41gbJAOf6AUlbuvAfKK+b5y8vxlytXfZKFp7z08A6yXzzyc/aL4PHCH3cc+7FUI5qJoLqfq9UOjqCqzwHPAWRlZTXYARS2bdvGuHHjWLRoEV6vl8zMTJ566il69OjBk08+yb333sv27dtp0aIFAPPmzeOCCy7ghRde4IYbbgBgyZIl9O/fn4kTJ3LXXXcxevRo5s+fT4sWLXC5XDzzzDOcddZZ0TxMY47fiVyQQ1yMDwmUK4WlfopKAxSVBSgsDZAy/0FSNi5gy3v3803f+ykqC1DsvEr85Ud8Likrp9gfoKSsnL9tGEaiVrkgZ79ICQmc7ZlKiT940S8NlFN1aJfV3l/jk7KK6ZTVr5Gy+jWKNYFTSqbUeqhP8ST3eV7nYnc2SVKKXwUQPFIOniQ4dQRc9Jdj+/vVIpqJIA/oWGk6A9gapVhOmKpy2WWX8etf/5pp06YBsHTpUrZv306PHj2YOnUqAwYM4N1332X06NEV2/Xu3Zs333yzIhFMmzaNvn37HvHdEydOZNSoUcyaNYubb76Z3NzcejsuY45SxxdzVaXEX87+ojIKSvwUFPspKC7jQLGfghI/l3/U/+hfx9kvUkoiV6S+R2GJn8LSAIWlfg6WBij1l1esWvVinLF+Khnrpx51MRYBn8eNN8FV8Z7odnFjyxe5sfBFBpUswEsJJeIlt9m5fNTuN1yUlEqi20Wix3X43eMiwe0iwS18XPYJ/df8lfbb5uAJFBNw+9iZcRGbsv7AW8lt8biEBLcLj1vwuILbeNwuPC7B4xKazV5E4jdfoW4vbn8JgoLHB4ES8DY/vjubGkQzEcwAbhWRacAgYJ+qHlUs1FjMnTuXhIQEbrnllop5p59+OgDr16/nwIEDTJw4kYceeuiIRNCpUyf279/P9u3bSUtL4+OPP+aSSy6pdh/nnXce69ati+hxGBNSDb/My8uV/cVl7D5Yyp7CUvYcLGNPYSn7isoY/dkgPDX8su5d9iqlgfJqdhQ0scqv42IS+TJhMG+0uIkWSQm0a+6jqddDk0S38/KQlOgiKdHDnMCnnL76cdK3zsYdKKbc7eNA1+HsP3c8X7dKx5fgxutcyGusU/tgLiz+HNw+vIFSBpzcmQEjhoTxh8qEfW1hayl4fLgDpbRNTaVt71PD+zsX74QzroOs62DaL4Pzrn4dsl+GA3VbYRzJ5qNTgfOBFBHJA+4HEgBU9e/ARwSbjq4j2Hz0urrY758/WMHKrfvr4qsq9GzfnPt/2qvWdZYvX84ZZ5xR7bKpU6dyzTXXcO6557JmzRp27NhBWlpaxfJRo0bx1ltv0a9fP/r374/X6632ez744AN69+59/AdiDBzTL/risgD5BSXsKCjh9Ckn4y4vObyw0sV8sHsqewpLKa+h4PZ5eZL7vVMZwtf4KKVEvHzT7Fxmd7yN65PTSfZ5aJ6UQHOfh2Sfh2RfAs28Hpp5PTT3JdBsTjbuxV+B24cvUMqQvicxZMSIMA62M+xOg7zgxdgVKKV5i9Y075QZ7l8LDu44fEE+1ovwiWx79euHP4+rVApQS7HY8YpYIlDVa0IsV+A3kdp/QzJt2jTeffddXC4Xl19+OW+99Ra/+c3hQ7/yyiu56qqrWL16Nddccw1ffvnlEdv/7ne/Y8KECaSmpvLiiy9W/Xpjjs38x9BNX1E8+yHWDvgzP+wrZvv+4or3HftL2LY/+Lmg2F+xWSpPVPll7uVr32Bmtv0NF7dMp03TRFo1SaR100RaNkmgVZPgdIsmCSR7Pbg+zIHFX1b8sh54cmcGjjg/vJgP5kfnYgxHXpCP9SJ8ItvWo5h4sriyUL/cI6VXr15Mnz79qPm5ubmsXbuWoUOHAlBaWkrXrl2PSATp6ekkJCTw6aefMmnSpKMSwaE6AmMqhPhVHyhXtu0vZvPuQvL2FJG3p5DffHE2CU7xjABJ30ymzzeT6eGUl3tcQlqyl7YtfHRLbcbgk9qQluwlLdlHarKXlGZeui78Gt/yw7/MzzutK+eNGBpezHX16zhGL8bRFHOJIFqGDBnCH/7wB55//nluuukmABYtWsTdd9/N+PHjuffeeyvW7dKlCxs3Htkj7AMPPMCOHTtwu931GrdppJxf9QWfTCDntD/x3c6DbNx1kO93FbJpdyF5ewopCxwupxGBWc3+j3tcr3Fm6QIStQS/y8eOjKHsO+d+vm7XkZSmXlyuEM+e/GdXdC7mJqIsEdQREeHdd99l3LhxPPLII/h8PjIzM5k3bx7PPnvkQ9OXXXYZ06ZNY9CgQRXzBg8eXN8hm2gLs6z+YImfdTsOsHbHAX42s19FpasAzZe/wgXLX+EsTSCL1+ncpgmntkvm4l7pdGrdhI6tk+jYqgntWvrwetzwwZfBik+PD0+glPZpabTv0T38mO1iHpNEqzZ+beCysrK06sA0q1at4tRTw6yJjzP2t2nAZt4JOS8Hf2GPeAJ/oJzvdh5k5Q/7WbOtgDXbCli9rYAte4sqNmnv3suDSdM4J7AQr/Orfl/mxXDxBFqndQz9NPm0X0Kztkf+oq98cTcxS0RyVDWrumV2R2BMfZuQBv6jW9/4NYGhTtt2j0vomtqU/p1bcfWAjnRvm0z3ts3o3LoJno+yYfF/Kn7Vt2mdAm3DfOLeftGbalgiMCbCysuVDTsPsmTTHpZu3svm5BcYtetZhkqw9U0RieQkncOiHnfyROcunJLenG5pzUj01NAV2Im2gjGmCksExpyIasr5i0oDLNm8h8Ub95C9Mfi+32mGmez10DsjjYyEdHw7ylC3F195Gef06sI5I84Nb5/2q97UMUsExpwIp/XO1vfH81rKWL7+bje5eXsrWux0T2vGJb3b0b9TK/p1aslJqc2CLXOmPQ2drrdf9aZBsERgzDHyB8pxPZSOKxAs5xegw7o3+P26NyglgSfO+Q8Du7TijE6tadGkhtHh7Fe9aUAsERgThm37ipm3Zgfz1uTzxbqd+Eqe4L6E1xnmzsFHCQG3Dz1lBInDHuKeOuwMzJj6YImgDuzatYsLL7wQCHZF7Xa7SU1NBWDGjBncfffdR3VNvXXrVh5//HFmzpwZzdBNDcrLlWVb9jFn1XZmr9rByh+C/Ve1a+HjJ33acW73Pvx4/SK8ucGnbN2BUkhqUac9QhpTX+I3EZzoQBWVtGnThqVLlwIwfvx4mjVrxl133YWqMnjw4Gq7pjYNiPO/4L/8JRbmJ/DJim18smIb2/eX4BLI6tyae4afwgUnp9GjbbPDbfVX7rTWOyYmxG8iCHOQixNRW9fU8+bN48CBA4waNaqi59LXXnsNESEnJ4c777yTAwcOkJKSwuTJk2nXrh3nn38+/fr1Iycnh/z8fF555RUefvhhli1bxlVXXcWECRMichyxzB8oJ3/GA7TduIB3nrydu4tH40tw8aMeqVzUM50hp6TRqmli9RtbOb+JEfGXCGp4mKcuh307pLauqSE4GtmKFSto3749Z599Nl988QWDBg3itttu4/333yc1NZU333yT++67j5deegmAxMREPv/8cyZNmsTIkSPJycmhdevWnHTSSdxxxx20adOmTo8hFqkq3+Tto9dL3UnQUg4NlH0ls7jSNwt1e5Ff1e3/gjENWfwlgrG58MkfYfVM8BdFZNi3cA0cOJCMjAwgeKfw/fff07JlS5YvX17RW2kgEKBdu3YV21x66aVAcGSzXr16VSzr2rUrmzdvtkRQi237inlnSR7Tc/LYkH+QDu5JPNFqOmcUfYEnUFzxvyBR+F8wJpriLxEkp4M3OTjcW4SGfTukpq6pD6k8AI3b7cbv96Oq9OrViwULFtS6jcvlOmJ7l8uF3++vdpt4VhYoZ86q7Uz9ejP/XptPucLALq25+byuDDutHS1mL4LFcyP+v2BMQxZ/iQDq7RH9mrqmLiwsrHGbk08+mfz8fBYsWMBZZ51FWVkZ3377Lb16RWechcZqy94iXv9qI//MzmPngRLatfBx6wXd+PkZGXRu0/TwitZdgzFxmgjqqZKvpq6pn3rqKbZs2VLtNomJiUyfPp3bb7+dffv24ff7GTdunCWC2jitfnTUy3y53cPkL79nzqrgBX3IKWn8YlAnftQjDXd1fe1bha8x1g11rIuHv41/xh24F7/MzMRh3Lb/V7RpmsjVAztyzcBOZLRqEu3wjGkQrBtqE5P0wTQkUFLxT/zT0n/xU9+/UPEiF1urH2PCVUM/t8Y0XPkFJTz8r1Vc4P8b7/kHUyLBSnP1JEHvK5Cxy6IcoTGNS8zcEahq6NGZ4kxjK/YLZcf+Yv7x+QZeX7iRUn85P+lzMue5uuJd/RV4fIi1+jHmuMREIvD5fOzatYs2bdpYMnCoKrt27cLn80U7lBO2+2Apz85bxysLNuIvV352egduHdKNLilNYdpea/VjzAmKicrisrIy8vLyKC4ujlJUDZPP5yMjI4OEhBq6Qm7gDpT4ef7zDbzw7w0UlQX4Wb8OjL2w+5HNP40xYYn5yuKEhAS6dOkS7TDMiajUCWBZk1SmLdrMpNnfsvNAKcNPS+e3F/WgW1pytKM0JibVmghEpH8Y31GmqlY7Z06M0wnglvfH86vtV7Eh/yADu7Tm+f86hX6dWkU7OmNiWqg7gvnAIoKDMNWkC5BZVwGZOFOlE8AO697gM94g0MSLa8x2q/Mxph6ESgSLVHVIbSuIyGd1GI+JM0X/vYTv3riDLvlzSZJSylw+XD1/ivviv4AlAWPqRa3PEYRKAuGuY0x15qzaztDn17B4mx+vlKFuLwlaittnTUCNqU8hK4tFJB1AVbeJSCpwLrBGVVdEOjgTm/ILShj/wQo+zP2B7mnNuKSrB1fa9dYE1JgoCVVZfDNwT/CjPAqMBlYAD4vIY6r6YuRDNLFCVXl78RYenLmSotIAd13UgzHnnUSi50eHV7KO34ypd6HuCG4FegFJwEagm3Nn0AqYC9SaCERkGDAJcAMvqOojVZZ3AqYALZ117lHVj47nQEzDtqOgmHvfXsac1TvI6tyKR37eh25pzaIdljGG0ImgTFULgUIRWa+q2wBUdY+I1Pokmoi4gWeAoUAesEhEZqjqykqr/RH4p6o+KyI9gY+wFkgxZ2buVv743nKKSgP8aURPrhuciau6LqGNMVERKhGUi0iCqpYBPzk0U0R8hO6wbiCwTlU3ONtMA0YClROBAs2dzy2ArccQu2ngDpT4uf/9Fby9OI++HVvy1yv62l2AMQ1QqERwOcGLNaqaV2l+G+C3IbbtAGyuNJ0HDKqyznhglojcBjQFflzdF4nIGGAMQKdOnULs1jQESzfvZey0JWzeXcjYC7tz25BueNzW2a0xDVGo5qObVLViIFwRaS4irYEiYHGI767u3r9qcdI1wGRVzQAuAV4VkaNiUtXnVDVLVbNSU1ND7NZETcE29OXhvD77a0Y9+yX+gPLmzWdxx9AelgSMacDC6mvIaT30AMEEcOhirkDXWjbLAzpWms7g6KKfG4BhAKq6wClySgFsVJFGqGTOIyRsXED5+kcZcsq9TBzVlxZNGmeHd8bEk3A7nbsL6KWqO4/huxcB3UWkC7AFuBr4RZV1NgEXApNF5FTAB+Qfwz5MQ+B0E+F1Jn/lmc2vNsyGJ7zwR8vpxjR04d6vrwcKj+WLnSKlW4FPgFUEWwetEJEHRORSZ7XfAjeJyDfAVGC0NrZ+sQ2fDp3FB+VnU0xicIYzUhg2UpgxjUK4dwT3Al+KyEKgoocwVb29to2cZwI+qjLvfyp9XgmcHXa0pkEJlCsTP1nD3+dv4R+tWuAt8oPHBzZSmDGNSriJ4B/AZ8AyoDxy4ZjG4kCJn7FTlzBn9Q6uPbMTPy5yIck2UpgxjVG4icCvqndGNBLTaGzdW8QNU7L5dnsBD/7sNH51Zmfg9cMrWDcRxjQq4SaCuU5b/g84smhod0SiMg3W8i37uG7yIopLA7w0egA/6mHNeY1p7MJNBIda+9xbaV6o5qMmxnz+bT7/77UcWjZJ5I3/HkT3tjZ0pDGxIFTvo+1U9QdVtQGB49zbOXn8/u1curdNZvJ1A2jb3BftkIwxdSTUHcFLTk+j84CPgf9UftLYxIcX/r2BCR+u4uxubfj7tWeQ7LOHxIyJJbUmAlUd7jztez5wGfC4iGwimBQ+VtVNkQ/RRIuq8tTstUyas5af9G7Hk1edTqLHuoowJtaErCNQ1WKcCz+A86TwcOB/RSRdVQdGNkQTDarKhA9X8eJ/vuPKrAwevrwPbus62piYFG5lcQVV/Q74P+D/RCSx7kMyUVWwDX3rOh5N/j0v5hzkurMz+dNPetr4AcbEsFCVxQUc7mTu0JVAnc+qqs2r3dA0WjrvUXTTAjr4n+aWH/2F3w87GRFLAsbEslB1BNY+MF44HccJwSz/K89sWDgbcqzjOGNiXdg1fyLSV0RudV59IhmUqX96+zd802ooRRos7VPrOM6YuBFWIhCRsQT7EEhzXq87o4qZGKCqPPrFPpbll+OVMtTjRazjOGPiRriVxTcAg1T1IICIPAosAJ6OVGCm/jwzdx1/n7+eD9uWId2vR6zjOGPiSriJQIBApekA1Q9FaRqZKV9+z+OzvuXyfh049YoZyKHWQdZxnDFxI9xE8DKwUETedaZ/BrwYmZBMfZmZu5XxH6xgaM+2PDaqjzURNSZOhZUIVPUJEZkHnEPwTuA6VV0SycBMZH25fid3vvkNWZ1b8fQ1/WxweWPi2LE8UPYd4He2ERHpr6qLIxOWiaSVW/cz5pUcMlOa8MJ/DcCX4I52SMaYKAorEYjIg8BogmMXH3rATIEhkQnLRMr2/cXcMGURzbweJl83kBZNrAM5Y+JduHcEVwInqWppJIMxkVVY6ufGKdnsLyrjrVsG075lUrRDMsY0AOEWDC8HWkYyEBNZ5eXKuGlLWbF1H0//oh8921vvIMaYoHDvCB4GlojIco4cqvLSiERl6tzEWWuYtXI79/+0J0NOsYfEjDGHhZsIpgCPAsuA8siFYyJhxjdbeXbeen4xqBOjB2dGOxxjTAMTbiLYqap/i2gkJiKWb9nH3dO/YUBmK8b/tJf1JGqMOUq4iSBHRB4GZnBk0ZA1H23Adh0o4eZXc2jVJJH/++UZNrqYMaZa4SaCfs77mZXmWfPRBixQrtw+bQk7D5Tw1i1nkZrsjXZIxpgGKtwniy+IdCCmjhRsg+nX8XzqH/li3V4eG9WHPhnW4MsYU7NaywpEZESoLwhnHVOP5j+GblxA06+e4MqsDK7M6hjtiIwxDVyoO4KJIrKF2nsafQiYWXchmePijDAGlUYYWz4bVtsIY8aY2oVKBNuBUP0Rr62jWMyJGJtL+Sf3Ubr8A3yUUO5JwnXqCLjoL9GOzBjTwIUas/j8eorDnKjkdJbll9NbSwm4vbhthDFjTJgi2p5QRIaJyBoRWSci99SwzpUislJEVojIG5GMJ5Z9tWEXP2zdzILWI3GPmQNnXGcjjBljwnIs3VAfExFxA88AQ4E8YJGIzFDVlZXW6Q7cC5ytqntEJC1S8cSyvYWl3PHmUnwt/sjMW84Br8dGGDPGhC2SdwQDgXWqusHptXQaMLLKOjcBz6jqHgBVtVrNY6Sq3PvOMnYeKOFvV/ejqTdiud0YE6NCJgIR6SwiKc7nM0XkLhG5LIzv7gBsrjSd58yrrAfQQ0S+EJGvRGRYDTGMEZFsEcnOz88PY9fx472lW/jX8m3cOfRkeme0iHY4xphGKNRzBH8CPgO+EpEJwFNACnC7iDwV4rura3KqVaY9QHfgfOAa4AUROerpJ1V9TlWzVDUrNTU1xG7jxw/7ivif91eQ1bkVY87rGu1wjDGNVKhyhGuAU4EmwCYgXVULRcQDLA2xbR5Q+WmmDGBrNet8paplwHcisoZgYlgUZvxxS1W5e3ou/oDy+BV9cdvA88aY4xSqaKhYVUtVdS+wXlULAVTVD4QarWwR0F1EuohIInA1wU7rKnsPuADAKX7qAWw4xmOIS699tZF/r93JfT85lcyUptEOxxjTiIW6I2gpIpcTLOZp7nzGma61QFpV/SJyK/AJ4AZeUtUVIvIAkK2qM5xlF4nISiAA/E5Vd53A8cSFzbsLefhfqzm3ewq/HNQp2uEYYxo5Ua1abF9pocjLtW2sqtfVeUQhZGVlaXZ2dn3vtsFQVf7rpa9ZvHEPn9xxHhmtmkQ7JGNMIyAiOaqaVd2yUE8W1/uF3tRuek4e/167kwdG9rIkYIypEzZSSSOyo6CYB2euZEBmK64d1Dna4RhjYoQlgkZk/IwVFPvLeeTnfXBZKyFjTB0J9RzBFc57l/oJx9Rk7uodfLRsG7cP6cZJqc2iHY4xJoaEuiO413l/O9KBmJoV7d5Cq3+OJKtNKWPOOyna4RhjYkyo5qO7RGQu0EVEqj4DgKpeGpmwTGWr3/wTfQOreLrDpyR6wundwxhjwheq+Wgi0B94Fbix6nJVnR+50KoXV81HK406dgSPjTpmjDk2tTUfrbVoyHmq+CtgsHPRXwzkqOr8aCSBeKO3f8MXSRdQRGJwhicJel8BY5dFNzBjTEwJt9VQWxFZAiwHVopIjoicFsG4DPDxRviuwIWPMvD4wEYdM8ZEQLiJ4DngTlXtrKqdgN8680yEFJcFmPDhKrr4CtEzrocbZ9uoY8aYiAh3FJOmqjr30ISqzhMR6+ksgv4xfwNb9hbhuul1XCe1Cc60UceMMREQbiLY4IxN8KozfS3wXWRCMlv2FvHs/HVc0judsw4lAWOMiZBwi4auB1KBd5xXCmD9EEXII/9ajSr84ZJTox2KMSYOhHVH4IwpfHuEYzHAkk17+OCbrdw+pJt1KmeMqRfW11ADovt/IPHVEZzctJCbf2RPEBtj6oclggZk83t/5tTSFTzdfhZNveFW3xhjzImxq01D4DxB3AlAoMfmf8L4f9oTxMaYehFWIhCRVOAmILPyNqp6fWTCijNjc9nw+jja/TCHJCkNPkF86gi46C/RjswYEwfCvSN4H/g3MJvg2MKmDh1MTCFnW4CfSxnq8SL2BLExph6FmwiaqOrvIxpJHHv5i+/oFtjDrl6/JPX8WyD7ZXuC2BhTb8JNBDNF5BJV/Sii0cShvYWl/OPzDZzZ/RGGXe10DGhPEBtj6lG4rYbGEkwGxSJS4Lz2RzKwuFCwjf3PXoSvZCe/vahHtKMxxsSpsBKBqiarqktVfc7nZFVtHungYl3h7IfJ2L+UJ9t+zCnp9uc0xkRH2M1HReRS4Dxncp6qzoxMSHHAaS7aBEDgnL0zYHwLay5qjImKsO4IROQRgsVDK53XWGeeOR5jcyk8+XKK1AacMcZEX7h3BJcAp6tqOYCITAGWAPdEKrCYlpzOip3lnEEZ6rbmosaY6DqWLiZaVvrcoq4DiSfb9xezJ38LC9uMRG6aYwPOGGOiKtw7goeBJSIyFxCCdQX3RiyqGPfsvPW86r+TedeeD62bWHNRY0xUhdsN9VQRmQcMIJgIfq+q2yIZWKzauXUjl2Rfj+u0h+nY2rqZNsZEX61FQyJyivPeH2gH5AGbgfbOPHOMvn/3frJkDeMS3412KMYYA4S+I7gTGAP8tZplCgypbWMRGQZMAtzAC6pabUsjERkFvAUMUNXsUEE3Sk6T0SwAgebLX4Hlr1iTUWNM1NWaCFR1jPNxuKoWV14mIr7athURN/AMMJTgncQiEZmhqiurrJdMcPSzhccYe+MyNpdVU24nM3+u9TBqjGlQwm019GWY8yobCKxT1Q2qWgpMA0ZWs96DwGNAcTXLYsbBUj9tdi7EK2Xg8YE1GTXGNBCh6gjSReQMIElE+olIf+d1PhCqprMDwfqEQ/KceZW/vx/QMdRTyiIyRkSyRSQ7Pz8/xG4bpo3vjCdV91LSsjvcONuajBpjGoxQdQQXA6OBDKByG8cC4A8htpVq5mnFQhEX8KTz/bVS1eeA5wCysrI0xOoNi1M30BNAIGnvt/D3c6xuwBjTYISqI5gCTBGRn6vq28f43XlAx0rTGcDWStPJwGnAPBEBSAdmiMilMVVhPDaXjVPvIG3LbKsbMMY0SOE+R/C2iPwE6AX4Ks1/oJbNFgHdRaQLsAW4GvhFpW33ASmHpp3nFO6KqSQAlJcrTbcuwGujjxljGqhwO537O3AVcBvBIp8rgM61baOqfuBW4BNgFfBPVV0hIg84PZnGhS0z/kwb3cOB5JOQG607CWNMwyOqoYvcRSRXVftUem8GvKOqF0U+xCNlZWVpdnYjuGlw6gaOYnUDxpgoEJEcVc2qblm4zRVdWCYAABAKSURBVEcPNe0sFJH2QBnQpS6Ci1ljc9l70kjratoY0+CF2+ncByLSEpgILCbY+uf5iEUVI+T7/+C1rqaNMQ1cyETgNPOco6p7gbdFZCbgcyp7TQ0KP32Y5v5d7PB1oe11r0L2y1Y3YIxpkEImAlUtF5G/Amc50yVANYXfBjhqGMq2Jd/ZcwPGmAYt3DqCWSLyc3Ea/JtajM3F32sURVjdgDGmcQi3juBOoCngF5Figk1IVVWbRyyyRqx07Tx8Wka5x4vL6gaMMQ1cuA+UJUc6kFih8x8lqXQn37s7kXnj65Az2eoGjDENWliJQETmqOqFoebFNadu4FDZWZfyTfCPc61uwBjT4IXqfdQnIq2BFBFpJSKtnVcm0L4+Amw0xubCaVdQKl4A1OoGjDGNRKg7gpuBcQQv+jkc7lF0P8FBZ8whyekUBcBXXoLflYDH6gaMMY1EqN5HJwGTROQ2VX26nmJqtEo2fIEPKO56Mc1atbW6AWNMoxBuZfHTIjIYyKy8jaq+EqG4GhenfqAlgECz9c44Ox5vNKMyxpiwhNv76KvA48A5wADnVW3nRXFpbC5bO46wfoWMMY1SuM8RZAE9NZyuSuNRcjqr98D5NuaAMaYRCvfJ4uUERxAz1diy+TtOL5jHytRLbMwBY0yjE+4dQQqwUkS+plI/Q6oaNwPM1CZ/5oP04SCJ6SmQ3htGPBF6I2OMaSDCTQTjIxlEo+VUEp8OwUriZVNg2RR7iMwY06iEVTSkqvOB74EE5/MiguMSxLexufxglcTGmEYu3FZDNwHTgX84szoA70UqqMbEt+XLioHpsUpiY0wjFG5l8W+Aswk+UYyqrgXSIhVUY3Fg1kO0DOxmd1IXqyQ2xjRa4dYRlKhq6aHhCETEQ3C4yvjk1A00AxBIKdpgg88YYxqtcO8I5ovIH4AkERkKvAV8ELmwGrixuehpoyi2wWeMMTEg3ERwD5APLCPYEd1HwB8jFVSDl5zOjtJEErWMgMvqBowxjVu4RUNJwEuq+jyAiLideYWRCqxBK9hG03Uf8j7nMfy6B3F/M8XqBowxjVa4iWAO8GPggDOdBMwCBkciqIau7LNHaBLYT3qbVvg69oWO9gCZMabxCjcR+FT1UBJAVQ+ISJMIxdRwOZXECQACZ+1+D8a3sEpiY0yjFm4dwUER6X9oQkTOAIoiE1ID5oxCVoKNQmaMiR3h3hGMA94Ska3OdDvgqsiE1IAlp3OgTGmqJfhdiTYKmTEmJoQ7MM0iETkFOJngcJWrVbUsopE1UP7vvwSgrNvFeFqkWSWxMabRC/eOAIKD0WQ62/QTkfgaoazKKGRJa53HKGwUMmNMI2cjlIVrbC67u15qHcwZY2JOREcoE5FhwCTADbygqo9UWX4ncCPgJ/jA2vWquvFY9lFvktNZt89FFmWo20YhM8bEjoiNUOY8dPYMMBzoCVwjIj2rrLYEyFLVPgR7N33sWPZRn8r2bqX7rtksbH4RcpN1MGeMiR2RHKFsILBOVTcAiMg0YCSwstL2cyut/xVwbZjx1LvtHzxAez1Ix7TWNgqZMSamRHKEsg7A5krTecCgWta/AfhXdQtEZAwwBqBTp07HEcoJcCqJMwAEMtZPhfFT7SEyY0zMOJYRylYDyc5rlTOvNlLdV1W7osi1BOshJtaw/+dUNUtVs1JTU8MJue6MzaWs58+tktgYE7PCbTV0JfA1cAVwJbBQREaF2CwP6FhpOgPYWnUlEfkxcB9wqaqWVF0edcnpbDroxksZ5W7radQYE3vCLRq6DxigqjsARCQVmE2wgrcmi4DuItIF2AJcDfyi8goi0o/g8JfDDn13g1OwjbTNH/Ox53yG3/gA5Ey2SmJjTEwJNxG4qlyodxHibkJV/SJyK/AJweajL6nqChF5AMhW1RkEi4KaEey+AmBTiAroelc8+2GaBgpo17YV0q6PVRIbY2KOhPNogIhMBPoAU51ZVwHLVPXuCMZWraysLM3Ozo78jpxK4qNYJbExphESkRxVrfZB4HAri39HsAinD9AXeC4aSaBeWU+jxpg4UWvRkIh0A9qq6heq+g7wjjP/PBE5SVXX10eQUZGczkFJIklLradRY0xMC3VH8BRQUM38QmdZ7CrYhmvV+7wdOJftV31oTxIbY2JWqMriTFXNrTpTVbNFJDMiETUU8x/D699P06bN6HDKQDhlYLQjMsaYiAiVCHy1LEuqy0AajEqVxC7gkpKPbDhKY0xMC1U0tEhEbqo6U0RuAHIiE1KUOZXEZa5gJXG5VRIbY2JcqDuCccC7IvJLDl/4s4BE4LJIBhY1yengTcZdXkopiSRaJbExJsbVmghUdTswWEQuAE5zZn+oqp9FPLIoOrB7K+/4L6TFOTcx0j/LKomNMTEt3DGL5wJzQ64YI95IHUvf9XdwUt+e0G5YtMMxxpiICndgmrihqqQv+RsDXGtIyYntFrLGGAPHNnh97JuQhvhLqOjsKPvF4MtaDBljYpjdEVQ2NpdVKRfb2APGmLhiiaASbdaWb/cKXikDj8/GHjDGxAVLBJWsWbeWc8v+w/cdLoUbZ1u3EsaYuGB1BJUc/PRhWnKQJqmtbIB6Y0zcsEQAFd1KnAEg4Fs6GZZOtkpiY0xcsKIhgLG57Op6qVUSG2PikiUCgOR0Nuxz4aUM9dgA9caY+GKJAPDv3Uq3XXNY2Pwi5MY5VklsjIkrVkcA7PhwAul6gA5WSWyMiUPxnQicSuL2AAKdNkyD8dOsktgYE1fiu2hobC6BXqMowiqJjTHxK74TQXI6W4o8eLWMgMsqiY0x8Sm+E0HBNlp9/xEfun6E3vCpVRIbY+JSXNcRFMx6iKaBAjq2a4OnQ1/oYJXExpj4E5+JwKkkTgYQOH3b2zD+baskNsbEpfgsGrJKYmOMqRCfiSA5nR/2l+HTUgKuRKskNsbEtfgsGgISti4EATl5ODRNsUpiY0zcir9E8GAaBEo49NtfVr0f/ODxRi0kY4yJpvgqGirYhnqTUcCvzqFb/YAxJs5FNBGIyDARWSMi60TknmqWe0XkTWf5QhHJjFgwD6bAX09GCncigEfKg/P9RVY/YIyJaxFLBCLiBp4BhgM9gWtEpGeV1W4A9qhqN+BJ4NGIBDMhDQJl1UUJfX9h9QPGmLgWyTqCgcA6Vd0AICLTgJHAykrrjATGO5+nA/8rIqKqWmdRjG9Ry0KBy56ts10ZY0xjFMmioQ7A5krTec68atdRVT+wD2hT9YtEZIyIZItIdn5+/rFFce274EniqMyS0AS6XXhs32WMMTEokolAqplX9Xoczjqo6nOqmqWqWampqccWRbch4E48ekdNU+Ha6cf2XcYYE4MimQjygI6VpjOArTWtIyIeoAWwu84j8ReDywMdBkBiMuAKDkBjjDEmonUEi4DuItIF2AJcDfyiyjozgF8DC4BRwGd1Wj9wyJ+s/yBjjKlJxBKBqvpF5FbgE8ANvKSqK0TkASBbVWcALwKvisg6gncCV0cqHmOMMdWL6JPFqvoR8FGVef9T6XMxcEUkYzDGGFO7+Hqy2BhjzFEsERhjTJyzRGCMMXHOEoExxsQ5iURrzUgSkXxg43FungLsrMNwGgM75vhgxxwfTuSYO6tqtU/kNrpEcCJEJFtVs6IdR32yY44PdszxIVLHbEVDxhgT5ywRGGNMnIu3RPBctAOIAjvm+GDHHB8icsxxVUdgjDHmaPF2R2CMMaYKSwTGGBPn4iYRiMgwEVkjIutE5J5oxxMJItJRROaKyCoRWSEiY535rUXkUxFZ67y3inasdUlE3CKyRERmOtNdRGShc7xvikhitGOsSyLSUkSmi8hq51yfFQfn+A7nf3q5iEwVEV+snWcReUlEdojI8krzqj2vEvQ353qWKyL9T2TfcZEIRMQNPAMMB3oC14hIz+hGFRF+4LeqeipwJvAb5zjvAeaoandgjjMdS8YCqypNPwo86RzvHuCGqEQVOZOAj1X1FKAvwWOP2XMsIh2A24EsVT2NYLf2VxN753kyMKzKvJrO63Cgu/MaA5zQ4OtxkQiAgcA6Vd2gqqXANGBklGOqc6r6g6oudj4XELxAdCB4rFOc1aYAP4tOhHVPRDKAnwAvONMCDAEOjUMaa8fbHDiP4FgeqGqpqu4lhs+xwwMkOSMZNgF+IMbOs6p+ztEjNNZ0XkcCr2jQV0BLEWl3vPuOl0TQAdhcaTrPmRezRCQT6AcsBNqq6g8QTBZAWvQiq3NPAXcD5c50G2Cvqvqd6Vg7112BfOBlpzjsBRFpSgyfY1XdAjwObCKYAPYBOcT2eT6kpvNap9e0eEkER41dD8Rsu1kRaQa8DYxT1f3RjidSRGQEsENVcyrPrmbVWDrXHqA/8Kyq9gMOEkPFQNVxysVHAl2A9kBTgkUjVcXSeQ6lTv/P4yUR5AEdK01nAFujFEtEiUgCwSTwuqq+48zefui20XmPlUGczwYuFZHvCRb3DSF4h9DSKUKA2DvXeUCeqi50pqcTTAyxeo4Bfgx8p6r5qloGvAMMJrbP8yE1ndc6vabFSyJYBHR3WhkkEqxomhHlmOqcUz7+IrBKVZ+otGgG8Gvn86+B9+s7tkhQ1XtVNUNVMwme089U9ZfAXGCUs1rMHC+Aqm4DNovIyc6sC4GVxOg5dmwCzhSRJs7/+KFjjtnzXElN53UG8F9O66EzgX2HipCOi6rGxQu4BPgWWA/cF+14InSM5xC8PcwFljqvSwiWm88B1jrvraMdawSO/XxgpvO5K/A1sA54C/BGO746PtbTgWznPL8HtIr1cwz8GVgNLAdeBbyxdp6BqQTrQMoI/uK/oabzSrBo6BnneraMYIuq4963dTFhjDFxLl6KhowxxtTAEoExxsQ5SwTGGBPnLBEYY0ycs0RgjDFxzhKBMcbEOUsEJm6JSBsRWeq8tonIlkrTX0Zgf6NFJF9EXqhlnSRn/6UiklLXMRhTHU/oVYyJTaq6i+DDWYjIeOCAqj4e4d2+qaq31hJTEXC6022GMfXC7giMqYaIHHDezxeR+SLyTxH5VkQeEZFfisjXIrJMRE5y1ksVkbdFZJHzOjuMffRyvmepM7hI90gflzHVsTsCY0LrC5xKsK/4DcALqjpQgiPA3QaMIzhYzJOq+h8R6QR84mxTm1uASar6utMHljtiR2BMLSwRGBPaInU69BKR9cAsZ/4y4ALn84+BnsE+0QBoLiLJGhwgqCYLgPucwXXeUdW1dR+6MaFZ0ZAxoZVU+lxeabqcwz+mXMBZqnq68+oQIgmgqm8AlwJFwCciMqSO4zYmLJYIjKkbs4CKSmAROT3UBiLSFdigqn8j2K1wn8iFZ0zNLBEYUzduB7KcSt+VBMv/Q7kKWC4iS4FTgFciGaAxNbFuqI2pJyIymmC/8TU2H6207vfOujsjHZcxdkdgTP0pAoaH80AZkECwDsKYiLM7AmOMiXN2R2CMMXHOEoExxsQ5SwTGGBPnLBEYY0yc+/+lOnXgzufr2QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.plot(campsol[:,0],campsol[:,5],label='CAMP')\n",
    "plt.plot(solTchem[:n,t_indx],solTchem[:n,C_indx],\"*\",label='TChem')\n",
    "plt.xlabel('Time [s]')\n",
    "plt.ylabel('Concentration of B [mol/m3]')\n",
    "plt.legend(loc='best')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
